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ABSTRACT 

We present the rest-frame optical (B, V, and R band) luminosity functions (LFs) of galaxies at 2 < z < 3.5, 
measured from a ^-selected sample constructed from the deep NIR MUSYC, the ultradeep FIRES, and the 
GOODS-CDFS. This sample is unique for its combination of area and range of luminosities. The faint-end 
slopes of the LFs at z > 2 are consistent with those at z ~ 0. The characteristic magnitudes are significantly 
brighter than the local values (e.g., ~ 1.2 mag in the R band), while the measured values for are typically 
~ 5 times smaller. The Z?-band luminosity density at z ~ 2.3 is similar to the local value, and in the R band 
it is ~ 2 times smaller than the local value. We present the LF of Distant Red Galaxies (DRGs), which we 
compare to that of non-DRGs. While DRGs and non-DRGs are characterized by similar LFs at the bright end, 
the faint-end slope of the non-DRG LF is much steeper than that of DRGs. The contribution of DRGs to the 
global densities down to the faintest probed luminosities is 14%-25% in number and 22%-33% in luminosity. 
From the derived rest-frame U—V colors and stellar population synthesis models, we estimate the mass-to-light 
ratios (M/L) of the different subsamples. The M/L ratios of DRGs are ~ 5 times higher (in the R and V bands) 
than those of non-DRGs. The global stellar mass density at 2 < z < 3.5 appears to be dominated by DRGs, 
whose contribution is of order ~60%-80% of the global value. Qualitatively similar results are obtained when 
the population is split by rest-frame U — V color instead of observed J—K color. 

Subject headings: galaxies: distances and redshifts — galaxies: evolution — galaxies: formation — galaxies: 
fundamental parameters — galaxies: high-redshift — galaxies: luminosity function, mass 
function — galaxies: stellar content — infrared: galaxies 



1. INTRODUCTION 

Understanding the formation mechanisms and evolution 
with cosmic time of galaxies is one of the major goals of ob- 
servational cosmology. In the current picture of structure for- 
mation, dark matter halos build up in a hierarchical fashion 
controlled by the nature of the dark matter, the power spec- 
trum of density fluctuations, and the parameters of the cosmo- 
logical model. The assembly of the stellar content of galaxies 
is governed by much more complicated physics, such as the 
mechanisms of star formation, gaseous dissipation, the feed- 
back of stellar and central supermassive black hole energetic 
output on the baryonic material of the galaxies, and mergers. 

The mean space density of galaxies per unit luminosity, or 
luminosity function (LF), is one of the most fundamental of 
all cosmological observables, and it is one of the most basic 
descriptors of a galaxy population. The shape of the LF re- 
tains the imprint of galaxy formation and evolution processes; 
the evolution of the LF as a function of cosmic time, galaxy 
type and environment provides insights into the physical pro- 
cesses that govern the assembly and the following evolution 
of galaxies. Therefore, the LF represents one of the funda- 
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mental observational tools to constrain the free parameters of 
theoretical models. 

The local (z ~ 0) LF has been very well determined from 
several wide-area, multi-wave band surveys with follow-up 
spectroscopy (Norberg et al. 2002; Blanton et al. 2001, 2003; 
Kochanek et al. 2001; Cole et al. 2001). At intermediate red- 
shifts (z ~ 0.75), spectroscopic surveys found a steepening 
of the faint-end LF with increasing redshift in the global LF, 
mainly due to the contribution by later type galaxies (Lilly 
et al. 1995; Ellis et al. 1996). From the COMBO-17 survey, 
Wolf et al. (2003) measured the rest-frame optical LF up to 
z ~ 1.2, finding that early-type galaxies show a decrease of 
a factor of ~ 10 in the characteristic density $* of the LF. 
The latest type galaxies show a brightening of ~ 1 mag in 
M* (the characteristic magnitude) and an increase of ~ 1.6 
in $* in their highest redshift bin in the blue band. Further 
progress in the measurement of the LF at z < 2 was obtained 
with the VIMOS VLT Deep Survey (VVDS; Le Fevre et al. 
2004) and the DEEP-2 Galaxy Redshift Survey (Davis et al. 
2003). From the VVDS data, Ilbert et al. (2005) measured 
the rest-frame optical LF from z = 0.05 to z = 2. From the 
same data set, Zucca et al. (2006) performed a similar anal- 
ysis for different spectral galaxy types, finding a significant 
steepening of the LF going from early to late types. Their 
results indicate a strong type-dependent evolution of the LF, 
and identify the latest spectral types as responsible for most 
of the evolution of the UV-optical LF out to z = 1 .5. 

Contrary to low-redshift studies, the selection of high- 
redshift (z > 2) galaxies still largely relies on their colors. 
One of the most efficient ways to select high-redshift galax- 
ies is the Lyman drop-out technique, which enabled Steidel 
and collaborators to build large samples of z ~ 3 star-forming 
galaxies (Steidel et al. 1996, 1999). Extensive studies of these 
optically (rest-frame ultraviolet) selected galaxies at z ~ 3 
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(Lyman Break Galaxies [LBGs]) and at z ~ 2 (BM/BX galax- 
ies; Steidel et al. 2004; Adelberger et al. 2004) have shown 
that they are typically characterized by low extinction, mod- 
est ages, stellar masses ~ 10 10 M©, and star formation rates 
of 10-100 Moyr" 1 (Steidel et al. 2003; Shapley et al. 2001; 
Reddy et al. 2005). Shapley et al. (2001) recovered the rest- 
frame V-band LF of LBGs at z ~ 3 from the rest-frame UV LF 
(Steidel et al. 1999; but see also Sawicki & Thompson 2006), 
finding that the LBG LF is characterized by a very steep faint 
end. 

LBGs dominate the UV luminosity density at z ~ 2 - 6, 
as well as possibly the global star formation rate density at 
these redshifts (Reddy et al. 2005). However, since the Ly- 
man break selection technique requires galaxies to be very 
bright in the rest-frame UV in order to be selected, it might 
miss galaxies that are heavily obscured by dust or whose light 
is dominated by evolved stellar populations. These objects 
can be selected in the near-infrared (NIR), which corresponds 
to the rest-frame optical out to z ~ 3. Using the NIR selection 
criterion J-K > 2.3 (also suggested by Saracco et al. 2001), 
Franx et al. (2003) and van Dokkum et al. (2003) discovered a 
new population of high-redshift galaxies (distant red galaxies 
[DRGs]) that would be largely missed by optically selected 
surveys. Follow-up studies have shown that DRGs constitute 
a heterogeneous population. They are mostly actively form- 
ing stars at z<~ 1.5-3.5 (van Dokkum et al. 2003; Ftirster 
Schreiber et al. 2004; Rubin et al. 2004; Knudsen et al. 2005; 
Reddy et al. 2005; Papovich et al. 2006). However, some 
show no signs of active star formation and appear to be pas- 
sively evolving (Labbe et al. 2005; Kriek et al. 2006; Reddy 
et al. 2006), while others seem to host powerful active galac- 
tic nuclei (van Dokkum et al. 2004; Papovich et al. 2006). 
Compared to LBGs, DRGs have systematically older ages 
and larger masses (Forster Schreiber et al. 2004), although 
some overlap between the two exists (Shapley et al. 2005; 
Reddy et al. 2005). Recently, van Dokkum et al. (2006) have 
demonstrated that in a mass-selected sample (> 10 11 M ) at 
2 < z < 3, DRGs make up 77% in mass, compared to only 
17% from LBGs (see also Papovich et al. 2006), implying that 
the rest-frame optical LF determined by Shapley et al. (2001) 
is incomplete. 

The global (i.e., including all galaxy types) rest-frame opti- 
cal LF at z > 2 can be studied by combining multiwavelength 
catalogs with photometric redshift information. Giallongo et 
al. (2005) studied the B-band LFs of red and blue galaxies. 
They find that the B-band number densities of red and blue 
galaxies have different evolution, with a strong decrease of 
the red population at z = 2- 3 compared to z — and a corre- 
sponding increase of the blue population, in broad agreement 
with the predictions from their hierarchical cold dark matter 
models. As all previous works at z > 2 are based on either 
very deep photometry but small total survey area (Poli et al. 
2003; Giallongo et al. 2005) or larger but still single field sur- 
veys (Gabasch et al. 2004), their results are strongly affected 
by field-to-field variations and by low number statistics, espe- 
cially at the bright end. Moreover, Gabasch et al. (2004) used 
an I-band-selected data set from the FORS Deep Field. The / 
band corresponds to the rest-frame UV at z ~ 2, which means 
that significant extrapolation is required. 

In this paper we take advantage of the deep NIR MUSYC 
survey to measure the rest-frame optical (B, V, and R band) 
LFs of galaxies at 2 < z < 3.5. Its unique combination of 
surveyed area and depth allows us to (1) minimize the effects 
of field-to-field variations, (2) better probe the bright end of 



the LF with good statistics, and (3) sample the LF down to 
luminosities ~ 0.9 mag fainter than the characteristic magni- 
tude. To constrain the faint-end slope of the LF and to in- 
crease the statistics, we also made use of the FIRES and the 
GOODS-CDFS surveys, by constructing a composite sample. 
The large number of galaxies in our composite sample also 
allows us to measure the LFs of several subsamples of galax- 
ies, such as DRGs and non-DRGs (defined based on their ob- 
served J-K color), and of intrinsically red and blue galaxies 
(defined based on their rest-frame U — V color). 

This paper is structured as follows. In § 2 we present 
the composite sample used to measure the LF of galaxies at 
2 < z < 3.5; in § 3 we describe the methods applied to measure 
the LF and discuss the uncertainties in the measured LF due to 
field-to-field variations and errors in the photometric redshift 
estimates; the results (of all galaxies and of the individual sub- 
samples considered in this work) are presented in § 4, while 
the estimates of the number and luminosity densities and the 
contribution of DRGs (red galaxies) to the global stellar mass 
density are given in § 5. Our results are summarized in § 6. 
We assume fi M = 0.3, = 0.7, and // = 70 km s" 1 Mpc -1 
throughout the paper. All magnitudes and colors are on the 
Vega system, unless identified as "AB". Throughout the pa- 
per, the J-K color is in the observed frame, while the U — V 
color refers to the rest frame. 

2. THE COMPOSITE SAMPLE 

The data set we have used to estimate the LF consists of 
a composite sample of galaxies built from three deep mul- 
tiwavelength surveys, all having high-quality optical to NIR 
photometry: the "ultradeep" Faint InfraRed Extragalactic Sur- 
vey (FIRES; Franx et al. 2003), the Great Observatories Ori- 
gins Deep Survey (GOODS; Giavalisco et al. 2004; Chandra 
Deep Field-South [CDF-S]), and the MUlti-wavelength Sur- 
vey by Yale-Chile (MUSYC; Gawiser et al. 2006; Quadri et 
al. 2007b). Photometric catalogs were created for all fields in 
the same way, following the procedures of Labbe et al. (2003). 

FIRES consists of two fields, namely, the Hubble Deep 
Field-South proper (HDF-S) and the field around MS 1054- 
03, a foreground cluster at z = 0.83. A complete descrip- 
tion of the FIRES observations, reduction procedures, and the 
construction of photometric catalogs is presented in detail in 
Labbe et al. (2003) and Forster Schreiber et al. (2006) for 
HDF-S and MS 1054-03, respectively. Briefly, the FIRES 
HDF-S and MS 1054-03 (hereafter FH and FMS, respec- 
tively) are K s band-limited multicolor source catalogs down 
to Kl ot = 24. 14 and = 23. 14, for a total of 833 and 
1858 sources over fields of 2.5' x 2.5' and 5.5' x 5.3', respec- 
tively. The FH and FMS catalogs have 90% completeness 
level at £ s tot = 23.8 and # s tot = 22.85, respectively. The final 
FH (FMS) catalogs used in the construction of the composite 
sample has 358 (1427) objects over an effective area of 4.74 
(21.2) arcmin 2 , with K%* < 23.14 (22.54), which for point 
sources corresponds to a 10 (8) o signal-to-noise ratio (S/N) 
in the custom isophotal aperture. 

From the GOODS/EIS observations of the CDF-S (data re- 
lease version 1 .0) a /f s band limited multicolor source catalog 
(hereafter CDFS) was constructed, described in S. Wuyts et al. 
(2007, in preparation). GOODS zero points were adopted for 
J and K s . The //-band zero point was obtained by matching 
the stellar locus on a J-K versus J-H color-color diagram 
to the stellar locus in FIRES HDF-S and MS 1054-03. The 
difference with the official GOODS //-band zero point varies 
across the field, but on average our //-band zero points are 
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TABLE 1 

The Composite Sample: Field Specifications 



Field 


Filter Coverage 


A 90 
(mag) 


K tot 

lim.cs 

(mag) 


Area 
(arcmin 2 ) 


N 




FIRES-HDFS 


U30oB45oV 60 6l814J s HK s 


23.80 


23.14 


4.74 


358 


68 


FIRES-MS 1054 


UBVV 6 o 6 I 814 J s HK s 


22.85 


22.54 


21.2 


1427 


297 


GOODS-CDFS 


1*435 V 60 6i775 Z850 JHK S 


21.94 


21.34 


65.6 


1588 


215 


MUSYC 


UBVRIzJHK s 


21.33 


21.09 


286.1 


5507 


116 



NOTE. — K$ is the K s band total magnitude 90% completeness limit, KfgJ, cs is the 7f s band total magnitude limit used to construct 
the composite sample, N is the total number of sources down to A^™ cs , and iV spec is the number of objects with spectroscopic redshift. 



~ 0.1 mag brighter. A total effective area of 65.6 arcmin 2 is 
well exposed in all bands. The final catalog contains 1588 
objects with #< ot < 21.34 in this area. At #< ot = 21.34 the 
median S/N in the K s isophotal aperture is ~ 12. 

MUSYC consists of optical and NIR imaging of four inde- 
pendent 30' x 30' fields with extensive spectroscopic follow- 
up (Gawiser et al. 2006). Deeper NIR JHK S imaging was 
obtained over four 10' x 10' subfiekls with the ISPI camera at 
the Cerro Tololo Inter- American Observatory (CTIO) Blanco 
4 m telescope. A complete description of the deep NIR 
MUSYC observations, reduction procedures, and the con- 
struction of photometric catalogs will be presented in Quadri 
et al. (2007b). The 5 u point-source limiting depths are 
J ~ 22.9, H ~ 21.8, and K s ~ 21.3. The optical UBVRIz data 
are described in Gawiser et al. (2006). The present work is re- 
stricted to three of the four deep fields: the two adjacent fields 
centered around HDF-S proper (hereafter MH1 and MH2) and 
the field centered around the quasar SDSS 1030+05 (M1030). 
The final MUSYC ^-selected catalog used in the construc- 
tion of the composite sample has 5507 objects over an effec- 
tive area of 286. 1 arcmin 2 , with Kl° l < 21 .09, which for point 
sources corresponds to a ~40 a S/N in the isophotal aperture. 

Table 1 summarizes the specifications of each field, includ- 
ing wave band coverage, Ks band total magnitude 90% com- 
pleteness limit (Kg°Q), effective area, the Ks band total mag- 
nitude limit used to construct the composite sample (K^ cs ), 
the number of objects, and the number of sources with spec- 
troscopic redshifts. 

Only a few percent of the sources in the considered catalogs 
have spectroscopic redshift measurements. Consequently, we 
must rely primarily on photometric redshift estimates. Photo- 
metric redshifts z p h ot for all galaxies are derived using an iden- 
tical code to that presented in Rudnick et al. (2001, 2003), but 
with a slightly modified template set. This code models the 
observed spectral energy distribution (SED) using nonnega- 
tive linear combinations of a set of eight galaxy templates. As 
in Rudnick et al. (2003), we use the E, Sbc, Scd, and Im SEDs 
from Coleman et al. (1980), the two least reddened starburst 
templates from Kinney et al. (1996), and a synthetic template 
corresponding to a 10 Myr old simple stellar population (SSP) 
with a Salpeter (1955) stellar initial mass function (IMF). We 
also added a 1 Gyr old SSP with a Salpeter IMF, generated 
with the Bruzual & Chariot (2003) evolutionary synthesis 
code. The empirical templates have been extended into the 
UV and the NIR using models. Comparing the photometric 
redshifts with 696 spectroscopic redshifts (63 at z ^ 1 -5) col- 
lected from the literature and from our own observations gives 
a scatter in Az/ (1 +z) of a = 0.06. Restricting the analysis 
to galaxies at z > 1.5 in the MUSYC fields gives a = 0.12, 
corresponding to Az m 0.4 at z = 2.5. Approximately 5% of 



galaxies in this sample are "catastrophic" outliers. A full dis- 
cussion of the quality of the photometric redshifts is given 
elsewhere (Quadri et al. 2007b). The effects of photometric 
redshift errors on the derived LFs are modeled in § 3.3. 

Rest-frame luminosities are computed from the observed 
SEDs and redshift information using the method extensively 
described in the Appendix of Rudnick et al. (2003). This 
method does not depend directly on template fits to the data 
but rather interpolates directly between the observed fluxes, 
using the best-fit templates as a guide. We computed rest- 
frame luminosities in the U, B, V, and R filters of Bessell 
(1990). For these filters we use M QjU = +5.66, M 0jB =+5.47, 
M Q] v = +4.82, and M 0:R = +4.28. In all cases where a spec- 
troscopic redshift is available we computed the rest-frame lu- 
minosities fixed at z S pec. 

Stars in all A" s -selected catalogs were identified by spec- 
troscopy, by fitting the object SEDs with stellar templates 
from Hauschildt et al. (1999) and/or inspecting their mor- 
phologies, as in Rudnick et al. (2003). On average, approxi- 
mately 10% of all the objects were classified as stars. 

We constructed a composite sample of high-redshift (2 < 
z < 3.5) galaxies to be used in the estimate of the LF in § 3. A 
large composite sample with a wide range of luminosities is 
required to sample both the faint and the bright end of the LF 
well; moreover, a large surveyed area is necessary to account 
for sample variance. The very deep FIRES allows us to con- 
strain the faint end of the LF, while the large area of MUSYC 
allows us to sample the bright end of the LF very well. The 
CDFS catalog bridges the two slightly overlapping regimes 
and improves the number statistics. The final composite sam- 
ple includes 442, 405, and 547 A" s -selected galaxies in the 
three targeted redshift intervals 2 < z < 2.5, 2.7 < z < 3.3, 
and 2.5 < z < 3.5, respectively, for a total of 989 galaxies with 
Kl ot < 23.14 at 2 < z < 3.5. Of these, -4% have spectroscopic 
redshifts. In Figure 1 we show the rest-frame B-band absolute 
magnitude versus the redshift for the composite sample in the 
studied redshift range 2 < z < 3.5. 

3. METHODOLOGY 

3.1. The 1/V max Method 

To estimate the observed LF in the case of a composite sam- 
ple, we have applied an extended version of the 1 /V max algo- 
rithm (Schmidt 1968) as defined in Avni&Bahcall (1980) so 
that several samples can be combined in one calculation. For 
a given redshift interval [21,22], we computed the galaxy num- 
ber density $(M) in each magnitude bin AM in the following 
way: 



1 N 1 
$(M) = ---£-, 
v ' AM r-( Vi 



(1) 
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FIG. 1 . — Rest-frame B-band magnitude vs. redshift for the composite sam- 
ple in the redshift range 2 < z < 3.5; the FH, FMS, CDFS, and MUSYC fields 
are plotted as red filled circles, blue open squares, orange open triangles, and 
green open circles, respectively. 



where N is the number of objects in the chosen bin and V; is: 



» rzupihJ) dV , 

j=\ J Zmm aZ 



(2) 



where wj is the area in units of steradians corresponding to 
the y'th field, n is the number of samples combined together, 
dV/dz is the comoving volume element per steradian, and 
Zup{i,j) is the minimum of Z2 and the maximum redshift at 
which the zth object could have been observed within the mag- 
nitude limit of the jth sample. The Poisson error in each 
magnitude bin was computed adopting the recipe of Gehrels 
(1986) valid also for small numbers. 

The 1 /V max estimator has the advantages of simplicity and 
no a priori assumption of a functional form for the luminosity 
distribution; it also yields a fully normalized solution. How- 
ever, it can be affected by the presence of clustering in the 
sample, leading to a poor estimate of the faint-end slope of the 
LF. Although field-to-field variation represents a significant 
source of uncertainty in deep surveys (since they are char- 
acterized by very small areas and hence small sampled vol- 
umes), the majority of published cosmological number densi- 
ties and related quantities do not properly account for sample 
variance in their quoted error budgets. Our composite sample 
is made of several independent fields with a large total effec- 
tive area of ~ 380 arcmin 2 (about a factor of 3 larger than the 
nominal area of the /^-selected CDFS-GOODS catalog used 
in Dahlen et al. 2005), which significantly reduces the uncer- 
tainties due to sample variance. Also, the large number of 
fields considered in this work with their large individual areas 
allows us to empirically measure the field-to-field variations 
from one field to the other in the estimate of the LF with the 
1 /Vmax method, especially at the bright end, and to properly 
account for it in the error budget. 

In order to quantify the uncertainties due to field-to-field 
variations in the determination of the LF, we proceeded as fol- 
lows. First, for each magnitude bin AM, we measured <£>■> for 
each individual jth field using equation (1). For each magni- 
tude bin with n 3, we estimated the contribution to the error 



budget of <£>(M) from sample variance using: 



(3) 



with n the number of individual fields used. For the magnitude 
bins with n ^ 2 (usually the brightest bin and the 3-4 faintest 
ones), we adopted the mean of the rms{&) with n 3. The 
final 1 a error associated to $(M) is then a = (cp oi + cr 2 v ) 1//2 , 
with erpoi the Poisson error in each magnitude bin. 

3.2. The Maximum Likelihood Method 

We also measured the observed LF using the STY method 
(Sandage et al. 1979), which is a parametric maximum likeli- 
hood estimator. The STY method has been shown to be unbi- 
ased with respect to density inhomogeneities (e.g., Efstathiou 
et al. 1988), it has well-defined asymptotic error properties 
(e.g. Kendall & Stuart 1961), and does not require binning 
of the data. The STY method assumes that $(M) has a uni- 
versal form, i.e., the number density of galaxies is separa- 
ble into a function of luminosity times a function of position: 
n(M,x) = $(M)p(x). Therefore, the shape of $(M) is deter- 
mined independently of its normalization. We have assumed 
that <£>(M) is described by a Schechter (1976) function, 

$(M) - (0.41nlO)$*[lO a4(M *- M)(1+Q) ] 

xexp[-10 a4 (">], (4) 

where a is the faint-end slope parameter, M* is the charac- 
teristic absolute magnitude at which the LF exhibits a rapid 
change in the slope, and $* is the normalization. 

The probability of seeing a galaxy of absolute magnitude 
M\ at redshift Zi in a magnitude-limited catalog is given by 



Pi 



$(Mi) 



(5) 



where Mf amt (z;) and Mbri g ht(zi) are the faintest and bright- 
est absolute magnitudes observable at the redshift Zi in 
a magnitude-limited sample. The likelihood A = I~Ii=i Pi 
(where the product extends over all galaxies in the sample) 
is maximized with respect to the parameters a and M* de- 
scribing the LF $(M). The best-fit solution is obtained by 
minimizing -2 In A. A simple and accurate method of esti- 
mating errors is to determine the ellipsoid of parameter values 
defined by 



In A = In A n 



(6) 



where x}j(S) lS the /3-point of the % 2 distribution with S de- 
grees of freedom. Parameter x^{S) is chosen in the standard 
way depending on the desired confidence level in the esti- 
mate (as described, e.g., by Avni 1976; Lampton et al. 1976): 
X«(2) = 2.3, 6.2, and 11.8 to estimate a-M* error contours 
with 68%, 95%, and 99% confidence level (1, 2, and 3 a, re- 
spectively). The value of $* is then obtained by imposing a 
normalization on the best-fit LF such that the total number of 
observed galaxies in the composite sample is reproduced. The 
1, 2, and 3 a errors on $* are estimated from the minimum 
and maximum values of $* allowed by the 1, 2, and 3 a con- 
fidence contours in the a-M* parameter space, respectively. 
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3.3. Uncertainties due to Photometric Redshift Errors 

Studies of high-redshift galaxies still largely rely on pho- 
tometric redshift estimates. It is therefore important to un- 
derstand how the photometric redshift uncertainties affect the 
derived LF and to quantify the systematic effects on the LF 
best-fit parameters. 

Chen et al. (2003) have shown that at lower redshifts (z < 1) 
the measurement of the LF is strongly affected by errors as- 
sociated with Zphot- Specifically, large redshift errors together 
with the steep slope at the bright end of the galaxy LF tend 
to flatten the observed LF and result in measured M* system- 
atically brighter than the intrinsic value, since there are more 
intrinsically faint galaxies scattered into the bright end of the 
LF than intrinsically bright galaxies scattered into the faint 
end. Using Monte Carlo simulations, Chen et al. (2003) ob- 
tained a best-fit M* that was 0.8 mag brighter than the intrinsic 
value in the redshift range 0.5 < z < 0.8. 

In order to quantify the systematic effect on the LF pa- 
rameters a and M* in our redshift range of interest (2 < 
z < 3.5), we performed a series of Monte Carlo simulations. 
The details of these simulations and the results are presented 
in Appendix A. Briefly, we generated several model cata- 
logs of galaxies of different brightness according to an in- 
put Schechter LF, extracted the redshifts of the objects from 
a probability distribution proportional to the comoving vol- 
ume per unit redshift (dV /dz), and obtained the final mock 
catalogs after applying a limit in the observed apparent mag- 
nitude. To simulate the errors in the redshifts, we assumed 
a redshift error function parametrized as a Gaussian distribu- 
tion function of 1 a width a' z (l+z), with a' z being the scatter 
in Az/(1 +z S pec), and we formed the observed redshift cat- 
alog by perturbing the input galaxy redshift within the red- 
shift error function. Finally, we determined the LF for the 
galaxies at z\ < z < Z2 using the 1 /V max and maximum like- 
lihood methods described in §§ 3.1 and 3.2, respectively. As 
shown in Appendix A, the systematic effects on the measured 
a and M* in the redshift interval 2 < z < 3.5 are negligible 
with respect to the other uncertainties in the LF estimate if 
the errors on the photometric redshifts are characterized by a 
scatter in Az/ (1 +z S pec) of a' z ~ 0.12, which is the appropri- 
ate value for the z > 2 sample considered in this work. This 
is not true at z < 1, where we find large systematic effects on 
both M* and a, consistent with Chen et al. (2003). As ex- 
plained in detail in Appendix A, the large systematic effects 
found at z < 1 arise from the strong redshift dependency of 
both dV /dz and dM /dz at low-z; at z > 2 these dependencies 
are much less steep, and this results in smaller systematic ef- 
fects on the measured LF. From the Monte Carlo simulations 
we also quantified that the effects of photometric redshift er- 
rors on the estimated luminosity density are typically a few 
percent (always < 6%). 

We conclude that the parameters of the LF and the lumi- 
nosity density estimates presented in this work are not signifi- 
cantly affected by the uncertainties in the photometric redshift 
estimates 8 . In order to include this contribution in the error 
budget, we conservatively assume a 10% error contribution to 
the luminosity density error budget due to uncertainties in the 
photometric redshift estimates. 

8 In Appendix A we also investigated the effects of non-Gaussian redshift 
error probability distributions. Systematic outliers in the photometric redshift 
distribution can potentially cause systematic errors in the LF measurements, 
although these are much smaller than the random uncertainties in the LF es- 
timates (if the outliers are randomly distributed). 
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FIG. 2. — Transmission curves of the J, H, and K s bands (gray shaded 
regions). The filled rectangles represent the observed wavelength ranges 
spanned by the central wavelengths of the rest-frame R band at 2 < z < 2.5 
(red), of the rest-frame V band at 2.7 < z < 3.3 (green), and of the rest-frame 
B band at 2 < z < 2.5 and 2.5 < z < 3.5 (blue; lower and upper rectangles, 
respectively). For the considered redshift intervals, the rest-frame R and V 
bands correspond approximately to the observed K s band. 

4. THE OBSERVED LUMINOSITY FUNCTIONS 

In this section we present the results of the measurement 
of the LF of galaxies at z > 2. We have measured the global 
LF in the rest-frame R and V band at redshift 2 < z < 2.5 
and 2.7 < z < 3.3, respectively. As shown in Figure 2, at 
these redshifts, the rest-frame R and V bands correspond ap- 
proximately to the observed K s band, which is the selec- 
tion band of the composite sample. We also measured the 
global LF in the rest-frame B band in the redshift interval 
2 < z < 2.5, to compare it with the rest-frame R -band LF, and 
at redshift 2.5 < z < 3.5, to compare it with previous studies. 
For each redshift interval and rest-frame band we also split 
the sample based on the observed J-K color (J-K > 2.3, 
DRGs; J-K < 2.3, non-DRGs) and the rest-frame U-V color 
{U-V > 0.25, red galaxies; U-V < 0.25, blue galaxies). In 
§ 4.1 we present the global LF of all galaxies, and in § 4.2 we 
present the LFs of the considered subsamples (DRGs, non- 
DRGs, red and blue galaxies) in the rest-frame R band. The 
results for the rest-frame V and B bands are shown in Ap- 
pendix B; in Appendix C we compare our results with those 
in the literature. 

4.1. Rest-Frame LF of All Galaxies 

Figure 3 shows the global rest-frame R- and B-band LFs 
for galaxies at 2 < z < 2.5, the rest-frame V-band LF at 
2.7 < z < 3.3, and the rest-frame B-band LF at 2.5 < z < 3.5. 
The large surveyed area of the composite sample allows the 
determination of the bright end of the optical LF at z > 2 with 
unprecedented accuracy, while the depth of FIRES allows us 
to constrain also the faint-end slope. This is particularly im- 
portant because of the well-known correlation between the 
two parameters a and M*. The best-fit parameters with their 
1, 2 and 3 a errors (from the maximum likelihood analysis) 
are listed in Table 2, together with the Schechter parameters 
of the local rest-frame B-band (from Blanton et al. 2003) and 
B-band (from Norberg et al. 2002) LFs. 

At redshift 2 < z < 2.5, the faint-end slope of the rest-frame 
B-band LF is slightly flatter than in the rest-frame B-band, al- 
though the difference is within the errors. In the two higher 
redshift bins, the faint-end slope of the rest-frame V-band LF 
is flatter (by ~ 0.3) than in the rest-frame B-band, although 
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Fig. 3. — Top left: Rest-frame B-band LF at 2 < z < 2.5 for all galaxies; black filled circles are the 1 /V max method estimates with 1 a error bars 
(Poisson errors only in black; field- to-field variation term included in gray); the red solid line is the LF estimated from the maximum likelihood 
method, and the red arrow represents the best-fit M*; the shaded region represents the 1 a uncertainties of the LF estimated from the maximum 
likelihood method; the brown dashed line and arrow represent the local B-band LF from SDSS (Blanton et al. 2003; R = r-0.12); the inset 
shows the best-fit value and the 1, 2, and 3 a confidence contour levels of the two parameters a and M*. Top right: Rest-frame V-band global 
LF at 2.7 < z < 3.3; the dark green solid line and the arrow represent the LF from the maximum likelihood method. Bottom left: Rest-frame 
6-band global LF at 2 < z < 2.5; the blue solid line and the arrow represent the LF from the maximum likelihood method; the black dotted line 
and arrow represent the local S-band LF from the 2dFGRS (Norberg et al. 2002; B = fej + 0.12). Bottom right: Rest-frame 6-band global LF at 
2.5 < z < 3.5; the cyan dashed line and arrow represent the rest-frame S-band LF of all galaxies at 2 < z < 2.5 for comparison; the black dotted 
line and arrow represent the local B-band LF from the 2dFGRS; in the inset, the 1 a confidence level of the estimated rest-frame fi-band LF at 
2 < z < 2.5 is also plotted. 



the difference is only at the 1 a level. Similarly, the faint- 
end slopes of the rest-frame B-band global LF in the low- and 
high-redshift bins are statistically identical. The characteris- 
tic magnitude Mg in the low-z interval is about 0.5 mag fainter 
with respect to the high-redshift one, although the difference 
is significant only at the ~1.5 a level. We therefore conclude 
that the rest-frame B-band global LFs in the low- and high- 
redshift bins are consistent with no evolution within their er- 



rors (< 2 a). 

In Figure 3 we have also plotted the local (z ~ 0.1) rest- 
frame /?-band (from Blanton et al. 2003) and B-band (from 
Norberg et al. 2002) LFs. The faint-end slope of the B-band 
LF at 2 < z < 2.5 is very similar to the faint-end slope of the 
local LF; the characteristic magnitude is instead significantly 
(> 3 er) brighter than the local value (by ~ 1.2 mag), and the 
characteristic density is a factor of ~ 4.8 smaller than the local 
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TABLE 2 

Best-fit Schechter Function Parameters for the Global LFs 







M* 






Redshift Range 


Rest-Frame Band 


(AB) 


a 


(10- 4 Mpc~ 3 mag" 1 ) 


2.0<z<2.5 
z = 


R 
R 


99 , 7 +0.20.0.35,0.45 
zz -° '-0.22.0.41,0.61 

-21.50±0.01 


, nl +0.21,0.35,0.49 
1 ' ul -0.20,0.34.0.47 

-1.05 ±0.01 


in ^c+2-65.4.50,5.92 
l"-"->_2.67,4.40,5.93 
51.1 =b 1.4 


2.7<z<3.3 
2.0<z<2.5 
2.5<z<3.5 
z = 


V 
B 
B 
B 


, 7 77 +0.20.0.33,0.45 
zz - ' '-0.24.0.43,0.63 
7 , Q , +0.22:0.39,0.51 
zl - y -0.26.0.49,0.73 
77 on+0.20;o.35,0.45 
-0.24.0.43,0.63 
-20.48 ± 0.07 


, . ,+0.23,0.40.0.57 
1 --0.25,0.41.0.57 
, .-+0.28,0.48:0.66 
1 - o -0.27,045.0.63 
, , n +0.25,0.42:0.58 
1 -^ u -0.24,041.0.56 
-1.24 ±0.03 


„ 7 +4.12,6.84,9.45 
,; >- u '-4.41 .7 .03 .9 .28 
1 A -.,+4.60.8.07,11.07 
'-4.64.7.53,9.84 
- fi +4.69:8.38, 11.49 
1J - ;,D -4.51,7.08,9.17 
62.9±3.1 



NOTE. — The quoted errors correspond to the 1, 2, and 3 a errors estimated from the maximum likelihood analysis as described 
in § 3.2. The z = values are the Schechter parameters of the local rest-frame fi-band (from Blanton et al. 2003) and B-band (from 
Norberg et al. 2002) LFs. We assumed Rab = r-0.12 and B A B = bj +0.12 to convert the local values into our photometric system. 



value. The rest-frame B-band LF at 2.5 < z < 3.5 is character- 
ized by a faint-end slope consistent with the local B-band LF; 
the characteristic magnitude is significantly brighter (> 3 a) 
than the local value by <~ 1.9 mag, while the characteristic 
density is a factor of ~ 4.6 smaller with respect to the local 
value. 

4.2. Rest-Frame LF ofDRGs, Non-DRGs, and Red and Blue 

Galaxies 

In this section we present the results of the LFs for different 
subsamples, by splitting the composite sample based on the 
observed J-K color (J-K> 2.3, DRGs; Franx et al. 2003) 
and on the rest-frame U—V color (by defining the red galaxies 
as those having U — V > 0.25, which is the median value of 
U-V of the composite sample at 2 < z < 2.5). In Figure 4, 
we show the rest-frame U—V color versus the observed J-K 
color for the composite sample at 2 < z < 3.5. 

In Figure 5 we show the rest-frame B-band LF at 2 < z < 2.5 
of DRGs versus non-DRGs and red versus blue galaxies, to- 
gether with the 1, 2, and 3 u contour levels in the a-M^ AB 
parameter space from the STY analysis. The LFs of the dif- 
ferent subsamples in the rest-frame V band at 2.7 < z < 3.3 
and in the rest-frame B band at 2 < z < 2.5 and 2.5 < z < 3.5 
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FIG. 4. — Rest-frame U-V color vs. observed J-K color for the composite 
sample at 2 < z < 3.5; the dashed vertical line represents the J-K > 2.3 
selection ofDRGs, while the horizontal dashed line shows the U-V > 0.25 
definition of red galaxies adopted in this paper. 



are shown in Appendix B in Figures. Bl, B2, and B3, respec- 
tively. In Table 3 the best-fit parameters and their 1, 2, and 3 a 
errors from the STY method are listed for all the considered 
rest-frame bands and redshift intervals. 

As shown in Figure 5, the rest-frame B-band LF at 2 < z < 
2.5 of DRGs is significantly (> 3 a) different from that of 
non-DRGs. The faint-end slope of the non-DRG LF is much 
steeper, indicating that the contribution of DRGs to the global 
luminosity and number density at faint luminosities is very 
small compared to that of non-DRGs. The bright end of the 
DRG LF is instead very similar to that of non-DRGs, with the 
two subsamples contributing equally to the global LF. Split- 
ting the composite sample based on the rest-frame U—V color, 
we find a qualitatively similar result, with the faint-end slope 
of the blue galaxy LF being much steeper than that of red 
galaxies (although the red galaxies clearly dominate the bright 
end of the LF). The difference between the LFs of DRGs (red 
galaxies) and non-DRGs (blue galaxies) is mainly driven by 
the different faint-end slopes. 

A similar result holds in the rest-frame V band at 2.7 < 
z < 3.3, although it is slightly less significant (at the 2-3 a 
level): the non-DRG (blue galaxy) LF is very similar to that of 
DRGs (red galaxies) at the bright end, while at the faint end, 
the LF of non-DRGs (blue galaxies) is steeper than that of 
DRGs (red galaxies). In the rest-frame B band, the differences 
between the LFs of DRGs/red galaxies and non-DRGs/blue 
galaxies become even less significant. Although DRGs/red 
galaxies are always characterized by LFs with flatter faint-end 
slopes, the significance of this result is only marginal (< 2 cr), 
especially in the higher redshift interval. 

Within our sample, there is marginal evidence for evolution 
with redshift: the rest-frame B-band non-DRG/blue galaxy 
LFs in the two targeted redshift bins are characterized by sim- 
ilar (within the errors) faint-end slopes, while the characteris- 
tic magnitude is brighter by <~ 0.5-0.6 mag in the higher red- 
shift bin. The LF of DRGs/red galaxies tends to get steeper 
from low to high redshifts and M* gets brighter by <~ 0.6 mag. 
However, because of the large uncertainties (especially for 
DRGs and red galaxies) on the measured Schechter param- 
eters, the differences in the rest-frame B band between the 
high- and the low-redshift bins are at most at the 2 a signifi- 
cance level. 

We note that the uncertainties on the estimated Schechter 
parameters mainly arise from the small number statistics at 
the very faint end, which is probed only by FIRES. Very deep 
(down to the deepest FIRES) NIR imaging over large spatially 
disjoint fields is required for further progress in our under- 
standing of the lowest luminosity galaxies at z > 2. 
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FIG. 5.— Rest-frame R-band LF at 2 < z < 2.5. Top: Rest-frame R-band LFs at 2 < z < 2.5 for DRGs/non-DRGs (left) and red/blue galaxies (right). The 
red solid line and filled circles represent the LF of DRGs and red galaxies estimated with the STY and 1 /V max methods, respectively; the blue dot-dashed line 
and filled circles represent the LF of non-DRGs and blue galaxies. The rest-frame R-band LF of all galaxies is also plotted (black dashed line). The shaded 
region (yellow for DRGs/red galaxies, cyan for non-DRGs/blue galaxies) represents the 1 a uncertainties of the LF measured with the STY method. The arrows 
represent the best-fit M*; error bars as in Fig. 3. Bottom: 1, 2, and 3 a contour levels from the STY method for DRGs and red galaxies (red dashed lines), for 
non-DRGs and blue galaxies (blue solid lines), and for all galaxies (black dot-dashed lines); the filled circles represent the best-fit Schechter values. The LF of 
non-DRGs (blue galaxies) is steeper than that of DRGs (red galaxies). 



4.3. Comparison with LBGs 

Shapley et al. (2001) computed the rest-frame optical (V 
band) LF of z ~ 3 LBGs using the distribution of optical if 
magnitudes (i.e., the rest-frame UV LF) and the distribution 
of M-K % colors as a function of Stf, magnitude. The rest- 
frame UV LF of LBGs was taken from Adelberger & Steidel 
(2000) with best-fit Schechter parameters a = -1.57 ±0.11, 
M* m = 24.54 ±0.14 mag, and $* = (1.5 ±0.4) x 10" 3 Mpc" 3 
in our adopted cosmology. Shapley et al. (2001) detected a 
correlation with 98% confidence between &-K s color and 2#, 
magnitude, such that fainter galaxies have redder M-K & col- 
ors. This trend was included in their LF analysis by using 
the relationship implied by the best-fit regression slope to the 
correlation, d(&-K s )/d& = 0.17 (the scatter around this 



regression is very large). The Schechter function was then 
fitted to the average LF values, obtaining best-fit Schechter 
parameters a = -1.85 ±0.15, M* AB = -22.99 ±0.25 mag, 

and = (6.2 ± 2.7) x 10" 4 Mpc" 3 . The overall shape of the 
rest-frame optical LF of LBGs is determined by the way in 
which the S#-K s distribution as a function of magnitude re- 
distributes magnitudes into K s magnitudes. Therefore, as 
a result of the detected positive correlation between Si, and 
S?-K s , the faint-end slope of the LBG rest-frame optical LF is 
steeper than that of the UV LF (Shapley et al. 2001). 

In Figure 6 we compare the rest-frame V-band LF of blue 
galaxies at 2.7 < z < 3.3 and the LBG LF from Shapley et 
al. (2001) in the same rest-frame band and redshift interval. 
The blue galaxy LF estimated with the 1/V ma x method ap- 
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TABLE 3 

Best-Fit Schechter Function Parameters for Subsamples 









M* 






Redshift Range 


Rest-Frame Band 


Sample 


(AB) 


a 


(10^* Mpc~ 3 mag" 1 ) 



2.0<z<2.5 



2.7<z<3.3 



2.0<z<2.5 



2.5<z<3.5 



J-K> 2.3 
J-K<23 
U-V > 0.25 
U-V <0.25 
J-K> 2.3 
J-K< 2.3 
>0.25 
C/-V<0.25 
/-A"> 2.3 
J-K<23 
U-V > 0.25 
U-V <0.25 
/-A"> 2.3 
7-^< 2.3 
J/-V >0.25 
C/-V<0.25 



-22.01 
-22.69; 

-22.09; 

-22.43 
-22.63 
-22.63 
-22.43 
-22.65 
-21.75 
-21.89; 
-21.63; 
-21.83 
-22.35 
-22.35 
-22.25 
-22.41 



26,0.43.0.57 
30,0.55.0.81 
27,0.43:0.57 
33,0.60.0.93 
20,0.35.0.47 
24,0.45.0.65 
32,0.51,0.69 
40,0.75.1.15 
28,0.47.0.63 
36,0.65.0.99 
24,0.41 ;o.55 
32,0.57.0.87 
28,0.45:0.61 
32,0.59.0.87 
26,0.43,0.57 
30,0.57.0.85 
42,0.67.0.87 
54,1.03,1.73 
26,0.43:0.57 
32,0.59.0.89 
30,0.49.0.65 
38,0.69.1.05 
32,0.51,0.67 
42,0.77.1.19 
32,0.51.0.67 
38,0.73.1.17 
26,0.43:0.57 
30,0.57.0.87 
32,0.51.0.67 
38,0.75.1.17 
26,0.43:0.57 
30,0.57.0.87 



+0.62.1.07.1.54 
-0.57.0.92.1.24 
_, , n +0.24,0.4l.0.57 
1 - 3U -0.24,0.40.0.55 
„ ,,+0.43,0.76.1.09 
U.JJ_ 44 _ 70 95 

+0.29,0.49.0.68 
-0.29,0.48.0.67 
+0.52,0.89.1.27 
-0.49,0.81.1.10 
+0.30,0.50:o.71 
-0.31,0.50.0.69 
+0.60,1.02.1.46 
-0.53,0.87.1.19 
+0.31,0.51.0.73 
-0.29,0.49.0.68 
+0.73,1.25.1.79 
-0.66,1.08.1.47 
+0.31,0.53:0.73 
-0.31,0.52.0.72 
+0.57,0.99.1.41 
-0.56,0.91.1.24 
+0.34,0.58:0.81 
-0.35,0.58.0.81 
+0.51,0.86.1.21 
-0.49,0.81.1.14 
+0.30,0.51.0.71 
-0.29,0.48.0.67 
+0.55,0.92.1.31 
-0.52,0.87.1.21 
+0.30,0.50.0.70 
-0.27,0.47.0.65 



0.62 



-1.39; 
-0.46; 
-1.21 
-0.11 
-1.26 
-0.54; 
-1.25 
-0.40; 
-1.33; 
-1.11 
-1.45 
-1.01 
-1.49 



+0.69,1.02.1.37 
1.09,1.90.2.58 
+2.72,4.59.6.30 
2.41,3.76.4.78 
+0.69,1.14.1.57 
1.14,2.19.3.31 
+2.97,4.97.7.32 
-2.42,3.64.4.47 
+1.10,1.81.2.56 
-1.71,3.09.4.34 
, , , .+4.02,7.12.9.88 
11 - 14 ^t.27,6.61.8.41 
fi 9 -+0.87,1. 54.2.24 
01 -0.92,2.08.3.40 
+4.50,7.72.10.55 
-4.16,6.52,8.36 
+1.29,2.18.3.18 
1.77,2.99.3.78 
+4.52,7.83.10.92 
-4.36,6.75.8.57 
+1.57,2.70.3.95 
2.33,4.35.6.20 
,+4.47.7.70,10.97 
-3.94.5.83,7.09 
+2.09.3.35.4.55 
-2.20,3.48.4.32 
,+4.53.7.95,11.13 



3.74 
6.32; 
7.51 
5.26; 
6.14" 



1 1.01 
4.12; 

11.08 
8.63 

8.19 



9.50; 
5.37" 



-3.92.6.04,7.57 
+2.05.3.25.4.49 
-2.34,3.77.4.73 



;.80" 



,+4.33.7.60, IOi 
-3.64.5.62,7.00 



NOTE. — The quoted errors correspond to the 1, 2, and 3 a errors estimated from the maximum likelihood analysis as described in § 3.2. 



pears consistent within the errors with the average LF values 
of LBGs (shown as stars in Figure 6). However, the best- 
fit Schechter parameters from the maximum likelihood analy- 
sis are only marginally consistent, with the faint-end slope of 
the LBG LF being significantly steeper than the one of blue 
galaxies, as shown in the inset of Figure 6. The same result 
is obtained if the rest-frame V-band LF of non-DRGs (rather 
than rest-frame blue galaxies) is compared to that of LBGs. 
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FIG. 6. — Rest-frame V-band LF of blue galaxies at redshift 2.7 < z < 
3.3 (blue dot-dashed line and filled circles) and z ~ 3 LBG LF (Shapley et 
al. 2001) (green dot-dashed and filled stars); the shaded area represents the 
1 £7 uncertainties of the blue galaxy LF measured with the STY method; the 
global LF measured with the STY method is also plotted as a black long- 
dashed line; arrows represent the best-fit M*. The inset shows the best-fit 
value and the 1,2, and 3 a confidence contour levels of the two parameters 
a and M* for blue galaxies (blue filled circle and solid lines) and the 68.3% 
confidence intervals for the LBG best-fit Schechter parameters from (Shapley 
et al. 2001) (green filled star and dashed line). 



4.4. Comparison with Previous Studies 

In Appendix C we compare our results with previously pub- 
lished LFs. Specifically, we have compared our rest-frame 
B-band LF with that derived by Poli et al. (2003), Giallongo 
et al. (2005), and Gabasch et al. (2004) in the redshift inter- 
vals 2 < z < 2.5 and 2.5 < z < 3.5, and our rest-frame /?-band 
LF with the rest-frame r-band LF derived by Gabasch et al. 
(2006) at 2 < z < 2.5. We also compared our rest-frame B- 
band LFs of red and blue galaxies at 2.5 < z < 3.5 with those 
measured by Giallongo et al. (2005). 



5.1 



5. DENSITIES 
Number Density and Field-to-Field Variance 



The estimates of the number density r\ obtained by inte- 
grating the best-fit Schechter function to the faintest observed 
rest-frame luminosity are listed in Table 4. For complete- 
ness, we also list rj Ml3SYC , calculated by integrating the best- 
fit Schechter LF to the rest-frame magnitude limits of the 
NIR MUSYC, and rf 2 , calculated by integrating the best-fit 
Schechter LF to 2 mag fainter than the faintest observed lumi- 
nosities. We find that the contribution of DRGs (red galaxies) 
to the total number density down to the faintest probed rest- 
frame luminosities is 13%-25% (18%-29%) depending on the 
redshift interval. By integrating the rest-frame B-band LF 
down to a fixed rest- frame magnitude limit (Mb.ab = -20), 
we find a hint of an increase of the contribution of blue galax- 
ies from the low redshift bin (62%) to the higher bin (74%), 
but the differences are not significant. If only the bright end 
of the LF is considered (integrating the LF down to the fixed 
NIR MUSYC limit, M b ,ab = -22.2), the increase of the con- 
tribution of the blue galaxy population becomes significant at 
the 2 a level, going from 42% in the low-z bin to 66% in the 
high-z bin. 

We determine the field-to-field variance in the density by 
fixing the parameters a and M* to the best-fit values measured 
using the composite sample, and estimating <P* for each y'th 
field separately by imposing a normalization on the LF such 
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TABLE 4 
Number Densities 









a log»? 


b logr? MUSYC 


c logr? +2 


Redshift Range 


Rest-Frame Band 


Sample 


(Mpcr 3 ) 


(Mpc- 3 ) 


(Mpc" 3 ) 



2.0<z<2.5 



2.7<z<3.3 



2.0<z<2.5 



2.5 <z<3.5 



All 

-«">2.3 
-K<23 
-V >0.25 
-V <0.25 
All 

-r>2.3 

-K<23 
-V >0.25 
-V <0.25 

All 
-K>23 
-K<23 
-V >0.25 
-V <0.25 

All 
-if > 2.3 
-£<2.3 
-V >0.25 
-V <0.25 



-2.61 
-3.48; 
-2.63 

-3.19; 
-2.72; 

-2.51 
-3.17; 
-2.63 
-3.25 
-2.61 
-2.53 
-3.34 
-2.60; 

-3.07; 

-2.71 
-2.43; 
-3.04 
-2.57 
-3.06; 
-2.56; 



+0.07.0.11,0.16 
-0.06.0.10,0.13 
+0.07.0.14,0.21 
-0.05.0.08,0.10 
+0.08.0.14,0.20 
-0.08.0.13,0.17 
+0.07:0.12,0.17 
-0.05.0.08,0.10 
+0.10:0.16,0.22 
-0.09.0.14,0.20 
+0.07:0.11,0.16 
-0.06.0.10,0.14 
+0.12:0.21,0.30 
0.10.0.16,0.20 
+0.08.0.13,0.18 
-0.07.0.12,0.16 
+0.1i:o.20,0.28 
-0.09.0.14,0.18 
+0.08:0.13,0.18 
0.07.0.12,0.16 
+0.08.0.13,0.19 
-0.08.0.12,0.17 
+0.17:0.30.0.42 
0.13.0.21,0.27 
+0.09.0.15,0.21 
-0.08.0.14,0.19 
+0.13:0.23.0.33 
0.11.0.17,0.22 
+0.10.0.16,0.23 
0.09.0.15,0.21 
+0.07.0.12,0.17 
-0.07.0.12,0.16 
+0.15:0.26,0.36 
-0.14.0.22.0.29 
+0.08:0.14,0.19 
0.08.0.13,0.18 
+0.16.0.27,0.38 
'-0.14.0.22.0.29 
+0.08:0.14,0.19 
-0.08.0.13,0.18 



-3.18; 

-3.64 
-3.37; 
-3.39; 
-3.59; 
-3.27; 
-3.65 
-3.49; 
-3.65 
-3.49; 
-3.38; 
-3.92; 
-3.53 
-3.63 
-3.71 
-3.46; 
-3.89; 
-3.65 
-3.90; 
-3.64 



+0.02.0.03,0.04 
-0.02.0.03,0.04 
+0.02.0.03,0.04 
-0.02.0.03,0.04 
+0.02:0.03,0.04 
-0.03.0.05,0.07 
+0.02.0.03.0.04 
-0.02.0.03,0.04 
+0.03.0.05,0.06 
-0.04.0.08,0.11 
+0.02:0.03,0.05 
-0.03.0.05,0.07 
+0.02:0.03,0.04 
-0.03.0.06,0.10 
+0.03.0.05,0.06 
-0.04.0.07,0.11 
+0.02.0.03,0.04 
-0.03.0.05,0.08 
+0.03.0.05.0.07 
-0.04.0.07,0.11 
+0.03.0.05,0.06 
-0.03.0.06,0.09 
+0.04.0.05,0.07 
0.05.0.11,0.17 
+0.04:0.06,0.07 
-0.05.0.08,0.12 
+0.03.0.04.0.06 
-0.04.0.08,0.12 
+0.04.0.07,0.09 
-0.06.0.10,0.16 
+0.03.0.05.0.06 
'-0.04.0.06,0.09 
+0.04:0.07,0.09 
-0.06.0.12,0.19 
+0.04.0.06,0.08 
-0.04.0.08,0.12 
+0.05.0.07,0.09 
-0.07.0.12,0.20 
+0.04:0.06,0.08 
-0.04.0.08,0.12 



-2.36; 
-3.48 
-2.24; 
-3.18 
-2.25 

-2.16; 

-3.05 
-2.22; 
-3.19; 
-2.17 
-2.15 
-3.19; 
-2.17 
-2.95 
-2.23 
-1.92; 
-2.68 
-2.03 
-2.74 
-2.00 



+0.17,0.29.0.42 
-0.13,0.20.0.27 
+0.09,0.19:0.31 
-0.06,0.09.0.11 
+0.25,0.42.0.59 
-0.18,0.29.0.39 
+0.09,0.17.0.28 
-0.06,0.09.0.11 
+0.31,0.52.0.73 
-0.21,0.34.0.46 
+0.19,0.33:0.47 
-0.15,0.24.0.32 
+0.26,0.50:o.76 
-0.17,0.25.0.31 
+0.24,0.41:0.59 
-0.19,0.30.0.40 
+0.21,0.40:o.64 
-0.13,0.19.0.24 
+0.24,0.41.0.59 
-0.20,0.32.0.41 
+0.25,0.42.0 60 
-0.17,0.28.0.37 
+0.42,0.80:i.23 
-0.23,0.33.0.40 
+0.30,0.51:0.73 
-0.20,0.33.0.43 
+0.29,0.58.0.88 
-0.18,0.26.0.32 
+0.35,0.59.0.83 
-0.22,0.37.0.49 
+0.23,0.39.0.54 
-0.18,0.30.0.40 
+0.42,0.76:1.10 
-0.32,0.47.0.58 
+0.27,0.46:0.65 
-0.22,0.35.0.47 
+0.44,0.78.1.14 
-0.31,0.45.0.56 
+0.26,0.45:0.64 
-0.22,0.35.0.47 



"The number densities are estimated by integrating the best-fit Schechter LF down to the faintest observed rest-frame luminosities: 
M R S AB = -19-6, = -20.3 and M B as ^ B = -19.5, M l ^ B = -20.0 for the low- and high-redshift intervals, respectively; they corre- 

spond to observed K s band magnitude of Kl ot ~ 23.14 at the lower limit of the targeted redshift intervals. The quoted errors correspond to 
the 1, 2, and 3 a errors estimated from the maximum likelihood analysis as described in § 3.2. 

b^MUSYC j s tng num b er density estimated by integrating the best-fit Schechter LF down to the rest-frame magnitude limits of the deep NIR 
MUSYC survey: M^| YC = -21.8, M^| YC = -22.4 and M^f YC = -21.7, -22.2 for the low- and high-redshift intervals, respectively 
(corresponding to observed A' ot ~ 21.09). 

c rf~ 2 is the number density estimated by integrating the best-fit Schechter LF down to 2 mag fainter than the faintest observed luminosities. 



that the total number of observed galaxies in each field is re- 
produced. In Table 5, the derived $* of DRGs and non-DRGs 
in each field are listed for the three targeted redshift intervals 
and compared to $* measured from the composite sample. 
The results in the redshift range 2.5 < z < 3.5 are plotted in 
Figure 7. We find an overdensity of DRGs in the M1030 field 
at all redshifts, with the excess (as compared to the character- 
istic density of the composite sample) varying from a factor 
of ~ 1 .2 in the lowest redshift bin up to a factor of ~ 1 .9 in the 
redshift interval 2.5 < z < 3.5. We also find an underdensity 
of DRGs (a factor of 0.82-0.86) in the GOODS-CDFS field, 
although only at z > 2.5. The value of for DRGs in M1030 
is a factor of ~ 1 .6-2.4 larger than that in the GOODS-CDFS 
field at z > 2.5, although they are similar at 2 < z < 2.5. 

These results are qualitatively consistent with van Dokkum 
et al. (2006), who showed that the GOODS-CDFS field is un- 
derdense in massive (M* > 10 11 M Q ) galaxies at 2 < z < 3, 
with a surface density that is about 60% of the mean and a fac- 
tor of 3 lower than that of their highest density field (M1030). 
However, our results seem to show systematically smaller un- 
derdensities for the GOODS-CDFS field compared to their 
work. In order to understand the origin of the smaller under- 
density of DRGs found for the GOODS-CDFS field in our 
work compared to that of massive galaxies in van Dokkum 
et al. (2006), we have estimated the surface density of DRGs 
in the redshift range 2 < z < 3 down to Kl° l = 21. We find 



that the surface density of DRGs in the GOODS-CDFS field 
is ~ 70% of the mean and a factor of ~ 2.2 ± 0.5 lower than 
that of the Ml 030 field, in good agreement with the values 




1 i. nDRG 

1 1 1 1 1 1 

FH FMS CDFS M1030 MH 1 MH2 
Field 

FIG. 7. — Values of for DRGs (red filled circles) and non-DRGs (blue 
filled triangles) in the redshift range 2.5 < z < 3.5 estimated in the rest-frame 
B band as described in § 5.1 for each individual field; the yellow (DRGs) and 
cyan (non-DRGs) shaded regions represent the values of <1>* from the analysis 
in § 3.2 on the composite sample. 
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TABLE 5 

Characteristic Density <I>* in Each Field 



Redshift Range 


Composite Sample 


MH2 


MH1 


Ml 030 


CDFS 


FMS 


FH 


DRGs 


2.0<z<2.5 
2.7<z<3.3 
2.5<z<3.5 


t 74 +0.69 
4.88!|;fe 


2.7 l +u " 
7.41+H 

8-54-1;! 


i 11+0.76 
1 75+155 


4.49^ 

8.49 + !;« 

9.37 + J:«« 


3 96 +u - 8J 
-Q.7.Q 

J - ZJ -0.81 
n nn+0.66 
J - yy -0.57 


2 70 +1 - 2S 
z -' u -0.88 

1 S4+ 074 


6 14+ 4 - !<4 

2 93 

7 7Q+3:84 
'• /y -2.69 
4 79+1.63 
-• /z -1.25 


non-DRGs 


2.0<z<2.5 
2.7<z<3.3 
2.5<z<3.5 


°- JZ -2.41 
11.14^1 


°- ZD -1.27 
1 1 40+3.80 
1 1 

7 rw+3.46 


8.63_! 09 

1 ? 44+ 2 - 44 

7 07+ 1 - 90 
'■"'-1.53 


7 90 +L0< > 

13 - 52+ 189 
14.15^1 


o 79+0.53 

y - D;, -1.03 
77+O.9I 
8 -''-0.83 


<: on+1.07 
3 - SU -0.91 
6.70+J-^ 

'- OJ -1.05 


12.65^-33 
18.68+"' 
12-68-111 



NOTE. — Units in 10 4 Mpc 3 mag the 1 a error of the characteristic density estimated for the individual field includes only the 
Poisson error (Gehrels 1986). The measured rest-frame R-, V-, and S-band LFs have been used in the redshift ranges 2.0 < z < 2.5, 
2.7 < z < 3.3, and 2.5 < z < 3.5, respectively. 



in van Dokkum et al. (2006). Therefore, the smaller under- 
densities of DRGs found for the GOODS-CDFS field in our 
work appear to arise mainly from the different targeted red- 
shift ranges. The approach adopted in this work to quantify 
field-to-field variance by comparing the $* of the individual 
fields might also mitigate field-to-field differences, especially 
at the bright end. 

We note that there are significant differences in the ob- 
served characteristic densities even within the MUS YC fields, 
although they have areas of ~ 100 arcmin 2 . For example, the 
observed $* of DRGs in the MH1 field is consistent with the 
one derived from the composite sample, but it is 0.61-0.76 
times the value in the Ml 030 field. These results demonstrate 
that densities inferred from individual ^100 arcmin 2 fields 
should be treated with caution. 

5.2. Luminosity Density 

In this section we present estimates of the luminosity den- 
sity. Because of the coupling between the two parameters a 
and M*, the luminosity density (obtained by integrating the 
LF over all magnitudes) is a robust way to characterize the 
contribution to the total LF from the different subpopulations 
and to characterize the evolution of the LF with redshift. 

The luminosity density p L is calculated using: 



PL 



L^L^dLv = T(2 + a)^*L*, 



(7) 



which assumes that the Schechter parametrization of the ob- 
served LF is a good approximation and valid also at lumi- 
nosities fainter than probed by our composite sample. Ta- 
ble 6 lists pl with the corresponding 1, 2, and 3 a errors 9 
for all of the considered samples. We also list p L ast , the lu- 
minosity density calculated to the faintest probed rest-frame 
luminosity, and p^ USYC , the luminosity density calculated to 
the rest-frame magnitude limits of the deep NIR MUSYC. 
While the difference between pl and p^ st is very small (neg- 
ligible for DRGs and red galaxies, and ^0.1 dex on average 
for non-DRGs and blue galaxies), the difference between p L 
and pMUSYC j s significant, especially for non-DRGs and blue 
galaxies (<~ 0.5 dex on average). 

9 The 1, 2, and 3 a errors of the luminosity densities were calculated by 
deriving the distribution of all the values of pl allowed within the 1,2, and 3 
a solutions, respectively, of the Schechter LF parameters from the maximum 
likelihood analysis. The contribution from the uncertainties in the photomet- 
ric redshift estimates derived in Appendix A was added in quadrature. 



In the top panel of Figure 8 we have plotted the total rest- 
frame B-band luminosity density versus the redshift, in- 
cluding a compilation of results from the literature. Only the 
results from the literature which are not significantly affected 
by field-to-field variations, or that have taken these into ac- 
count, are plotted. Our measurement of the total rest-frame 
B-band luminosity density is the only one at z > 2 that is 
not significantly affected by field-to-field variance. From Fig- 
ure 8, there is an indication of a possible increase of the total 
luminosity density in the highest redshift bin, significant at 
the ~ 2.2 a level. The measurement at 2 < z < 2.5 is con- 
sistent with the one at z ~ 1.9 from Dahlen et al. (2005). In 
Figure8 we have also plotted the computed B-band rest-frame 
luminosity density as a function of z predicted from large- 
scale ACDM hydrodynamical simulations from Nagamine et 
al. (2000, 2001) and from a semianalytical model taken from 
Dahlen et al. (2005). While the predicted luminosity densities 
match the measurements at z < 1.2 well, they clearly over- 
predict them at larger redshifts. Only the prediction at z ~ 3 
from Nagamine et al. (2000) is consistent with our measure- 
ment, although their model still overpredicts significantly the 
luminosity densities in the range 1 .2 < z < 2.5. 

The bottom panel of Figure 8 shows the total rest-frame R- 
band luminosity density p L versus the redshift, including a 
compilation of results from the literature as in the top panel. 
As for the rest-frame B band, our measurement of the rest- 
frame B-band p^ is the first one at z > 2 for which sample 
variance does not significantly contribute to the error budget. 
Our point at z > 2 is consistent with the trend observed at 
z < 2 of decreasing luminosity densities with increasing red- 
shifts, although the point at z ~ 1.9 from the LF analysis of the 
GOODS-CDFS survey by Dahlen et al. (2005) is only ~ 43% 
of our measurement at z ~ 2.2. From Table 5 we see that in 
the redshift range 2 < z < 2.5, the CDFS field is underdense 
in non-DRGs. The value of <f>* for all galaxies in the CDFS 
field is 6.67+[] f 4 x 10" 4 Mpc" 4 mag" 1 (estimated as in § 5.1), 



which is a factor of 0.63 



+0.30 
-0.18 



the value of from the com- 



posite sample (see Table 2). Therefore, the lower value of the 
rest-frame B-band p L from Dahlen et al. (2005) at 1 .8 < z < 2 
could be due to an underdensity of galaxies at z ~ 2. In Fig- 
ure 8 we have also plotted the computed B-band rest-frame 
luminosity density as function of z predicted from large-scale 
ACDM hydrodynamical simulations from Nagamine et al. 
(2000). As for the rest-frame B band, the predictions match 
the observations well enough at z < 1, but at larger redshifts 
they significantly overpredict them. 
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TABLE 6 
Luminosity Densities 



Redshift Range Rest-Frame Band Sample logpL logp^ s ' a logp^ USYC 

7 fl <r <r 9 S )? All 9 , .,.+0.06,0.08.0.12 9 , .,,+0.05,0.06.0.07 ,,,-,+0.05.0.06,0.07 

z - u — z — K An ZD - ' -0.05,0.06.0.07 ZD - ' -0.05,0.06.0.07 zo --' ; '-0.05.0.06,0.07 

9q 9 , . ,+0.05,0.06.0.08 9f - , „+0.05,0.06;o.08 ? , +0.05:0.07,0.09 

R > Z J ZD - 1 -0.05,0.06.0.07 ZD - 1 -0.05,0.06.0.07 Z0 - 1 -0.05.0.07,0.09 

,_ „ , 9 , ,, , .+0.10,0.20:o.43 9 , , 7 +0.05,0.06;o.07 9 , ,,+0.06:0.08,0.10 

J R - Z - J ZD - o4 -0.07,0.08.0.10 Z0 - : > '-0.05,0.06,0.07 ZO -' U -0.06.0.09,0.11 

r ,_, />fl9 c 9 , .,+0.05,0.06.0.07 9 , .,+0.05,0.06:0.07 ? , , 7 +0.05:o.06,0.07 

(7 V s U.Z3 ZD - 4: '-0.05,0.06.0.07 ZD - -0.05,0.06.0.07 M i -0.05.0.06,0.07 

j, n «, 9 , ,,+0.18,0.48:2.47 9 , .,+0.05,0.06:0.08 ,,,,+0.07,0.10,0.13 

U V < U.Z3 ZD - :,U -0.09,0.12.0.14 ZD - 4u -0.05,0.06.0.07 ZO - U:, -0.08.0.12,0.16 

„ ... 9 , Q7 +0.08,0.15;o.26 9 , Q1 +0.05,0.06;o.07 ,. +0.05:0.07,0.08 

Z - / - z - J - J V AU Zt, - y -0.06,0.07.0.09 Z "- y -0.05,0.06.0.07 Z "- t,1 -0.06.0.08,0.10 

9q 9 , ..+0.07,0.12:0.24 9 , .,+0.05,0.06:0.08 9fi +0.06:o.07,0.09 

R > Z -J ZD - 44 -0.05,0.06.0.07 ZD - -0.05,0.06.0.07 ZO - ZO -0.06.0.09,0.13 

, „ , 9 - „ 9 +0. 12,0.28:0.66 9 , 7 o+0.05,0.06:o.07 -,+0.06,0.08,0.11 

7 A. S Z -J Zt) - 8Z -0.07,0.10.0.12 Zt> - '•'-0.05,0.06.0.07 ZO - ; ' O -0.07,0.10,0.13 

11 -V ■> riic 9<; An +0.05,0.08.0.14 9 , ,,+0.05,0.06:0.07 9A 9 ^+0.06.0.07,0.09 

(7 V s U.Z3 ZD - 4u -0.05,0.06.0.07 ZD - Jy -0.05,0.06.0.07 

„ 9 , „.+0.13,0.3i;o.85 9 , _ .+0.05,0.06;o.07 

ZD -° 4 -0.08,0.11.0.13 za /4 -0.05,0.06.0.07 ^ U - Ju -0.07.0.10,0.13 

, n , „ ... 9 , ,.+0.10,0.20:o.38 9 , -.+0.05,0.06:0.07 ,,.,+0.06,0.08,0.11 

Z -USZS Z -J d All Zt) - t,1 -0.07,0.09.0.11 lbi -0.05,0.06.0.07 ZO - I;, -0.06.0.09,0.13 

9 , 9 , Qn +0.10,0.27:2.95 9< - SQ +0.06,0.09:0.13 ,,- +0.08:0.11,0.15 

J K > Z.i z:> - yu -0.06,0.07.0.08 Z:> - 5y -0.05,0.06.0.07 Z -'- -0.09.0.15,0.22 

, v ^ 9 , ,,-,+0.15,0.35,1.15 9 , .,+0.05,0.07:0.08 ? - QQ +0.07:0.10,0.13 

J K S Z -J Zt>.r>jS_ og j j 13 zt, - 4;, -0.05,0.06.0.07 Z;, - yy -0.07.0.12,0.16 

rj-i/^ms ,,,,+0.07,0.14,0.31 9 , , 7 +0.05,0.07:0.10 Q9 +0.06:o.08,0.11 

(7 V £ U.Z3 ZD - 1 °-0.05,0.06.0.07 Z0 ' 1 -0.05,0.06.0.07 Z:, - yZ -0.07.0.11,0.16 

r7 „ n ,c 9 , .,+0.22,0.68:2.38 9 , 9Q +0.05,0.07:0.09 9 - 7Q +0.08:0.12,0.16 

U V < U.ZD ZD - 41 -0.10,0.14.0.16 ZD - zy -0.05,0.06.0.07 /y -0.09.0.14,0.21 

9 . „ ... Qn +0.16,0.37:i.OO 9 , 7 ,+0.05,0.06:o.08 ,, 9 „+0.06:o.08,0.09 

ZJ<Z - J C A11 ZD - yU -0.10,0.13.0.16 Z0 - /;, -0.05,0.06.0.07 ZO ' ZB -0.06.0.09,0.12 

T v -.90 9 , ,,+0.21,0.81,°° 9 , 9 ,+0.07,0.10:0.14 9 c o fi +0.07:0.10,0.12 

ZD - -0.10,0.13.0.14 ZD - ZJ -0.06,0.07.0.09 ^ -0.09.0.15,0.22 

, 9 , 9 , .,,+0.24,0.74,°° 9 , -„+0.05,0.07:o.09 ,, n „+0.07:o.09,0.11 

J HviJ ZD. '0-0.12,0.16.0.19 ztoy -0.05,0.06.0.07 ZO - U °-0.07.0.11,0.15 

r ,_ v>n9 c 9 , 97 +0.20,0.67,«o 9 , 99 +0.07,0.10:0.14 ,,,,1+0.07,0.10,0.12 

U V s_ u.ZJ ZD - Z/ -0.09,0.12,0.13 ZD - zz -0.06,0.07,0.09 -0.09,0.15,0.23 

t/-V<0.25 iu - /;7 -0.13,0.17,0.20 iu - uu -0.05,0.06,0.07 ^ u ' u7 -0.07 ,0.10,0.14 



NOTE. — The luminosity densities are in units of erg s Hz Mpc . The quoted errors correspond to the 1, 2, and 3 a errors estimated 
from the maximum likelihood analysis as described in § 3.2 and include the contribution from photometric redshift uncertainties. 

a PL St is the luminosity density estimated by integrating the LF down to the faintest observed luminosities: M^^ B = -19.6, My S ^ B = -20.3 
and Mg S ^ B = -19.5, M B a ^ B = -20.0 for the low- and high-redshift intervals, respectively; they correspond to observed K s band magnitude 
of Kl° l ~ 23. 14 at the lower limit of the targeted redshift intervals. 

bpMUSYC - s luminosity density estimated by integrating the LF down to the rest-frame magnitude limits of the deep NIR MUSYC sur- 
vey: M^| YC = -21.8, M^USYC = _ 22 .4 an d M B yI ^§ YC = -21 .7, -22.2 for the low- and high-redshift intervals, respectively (corresponding 
to observed A: s tot ^ 21.09). ' 



We can also compare our estimated luminosity densities 
in the rest-frame B and V bands with those from Rudnick 
et al. (2006), who presented the evolution of the rest-frame 
optical luminosity and stellar mass densities at z < 3. The 
luminosity density in Rudnick et al. (2006) was computed 
by simply adding up the luminosities of all galaxies in the 
targeted redshift bins with rest-frame V-band luminosities 
> 3 x 10 10 L©. If we integrate our measured LF down to 
the same limit adopted in Rudnick et al. (2006), we obtain 
log pl = 26.85 ± 0.06 and 26.66 ±0.06 for the V and B band, 
respectively, at z ~ 3, and 26.48 ±0.06 for the B band at 
z ~ 2.2 (units in erg s _1 Hz -1 Mpc -3 ), in excellent agree- 
ment with their estimates (26.79 ± 0.05, 26.66 ± 0.05, and 
26.37 ± 0.04, respectively). 

In § 4.2 we showed that the contribution of DRGs and red 
galaxies to the global LF is comparable to (or larger than) 
that of non-DRGs and blue galaxies at the bright end, but 
it becomes significant smaller at the faint end, where non- 
DRGs and blue galaxies dominate the global LF. In § 5.1, 
the contribution of DRGs to the global number density has 
been shown to be 13%-25% down to the faintest probed rest- 
frame luminosities. The contribution of DRGs to the global 
luminosity density is 19%-29% depending on the considered 
rest-frame band and redshift interval (see Table 7). Their con- 
tribution increases up to 31%-44% if we cut the composite 
sample to the rest-frame absolute brightness limit of MUSYC 



(~ 0.45L*, 0.7L*, 0.8L* of the global LF in the rest-frame 
R, V, and B band, respectively), which reflects the increas- 
ing importance of DRGs at the bright end. A similar result 
holds if we consider the red galaxy subsample. Their contri- 
bution to the global luminosity density is 29%-52% down to 
the faintest observed luminosities and increases up to 36%- 
69% if we limit the analysis to the MUSYC brightness limits. 
From Table 7 we conclude that the total luminosity density 
is dominated by non-DRGs/blue galaxies, especially in the 
bluer rest-frame optical bands, although DRGs/red galaxies 
contribute about 50% at the bright end. 

5.3. Stellar Mass Density Estimates 

Although non-DRGs and blue galaxies represent the major 
contribution to the total luminosity and number densities in 
the rest-frame optical bands, it has been shown that DRGs 
usually have larger mass-to-light ratios than non-DRGs (e.g., 
Forster Schreiber et al. 2004; Labbe et al. 2005; van Dokkum 
et al. 2004), as is generally true for red versus blue galaxies 
(e.g., Bell & de Jong 2001; Kauffmann et al. 2003). It is 
therefore interesting to quantify the contribution of DRGs (red 
galaxies) to the total stellar mass density. 

Following the method described in Rudnick et al. (2003), 
we estimated the stellar mass density from the measured 
global luminosity density modulo the mass-to-light ratio 
M/L. For each subsample we have measured the median rest- 
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Fig. 8. — Top: Total rest-frame B-band luminosity density pl plotted vs. the redshift. Blue filled stars are the values from this work (the error 
bars include also the contribution from photometric redshift uncertainties); thick error bars represent statistical errors only, while thin error bars 
also include the sample variance contribution to the error budget. The z < 2 points are taken from the 2dFGRS (Norberg et al. 2002; asterisks), 
Autofib Redshift Survey (Ellis et al. 1996; filled squares), COMBO- 17 (Wolf et al. 2003; open squares), VIMOS-VLT Deep Survey ((llbert et 
al. 2005); open diamonds), and GOODS-CDFS (Dahlen et al. 2005; B- and A"-selected catalog as open and filled circles, respectively). The 
long- and short-dashed lines correspond to the computed B-band rest-frame luminosity densities as a function of redshift using large-scale 
ACDM hydrodynamical simulations from Nagamine et al. (2000) and Nagamine et al. (2001), respectively; the dot-dashed line is a prediction 
from a semianalytical model taken from Fig. 13 of Dahlen et al. (2005). Bottom: Total rest-frame B-band luminosity density pl plotted vs. the 
redshift. The red filled star represents the value from this work. The z < 2 points are taken from the CNOC1 Redshift Survey (Lin et al. 1997; 
filled squares), Century Survey (Brown et al. 2001; filled triangle), SDSS (Blanton et al. 2001; asterisk), COMBO-17 (Wolf et al. 2003; open 
squares), Las Campanas Infrared Survey (Chen et al. 2003; open circles), GOODS-CDFS (Dahlen et al. 2005; filled circles), and VIMOS-VLT 
Deep Survey (llbert et al. 2005; open diamonds). The long-dashed line corresponds to the computed B-band rest-frame luminosity density as a 
function of redshift using large-scale ACDM hydrodynamical simulations from Nagamine et al. (2000). 



frame U — V color 10 , estimated the corresponding M/L ratio 
from the relation between U—V color and M/L ratio obtained 
from stellar population synthesis models, and multiplied the 
estimated M/L ratio by the measured luminosity density to 
obtain the stellar mass density. 

10 Rather than the median U — V color of the sample, Rudnick et al. (2003) 
computed the global U — V color from the relation: U — V = -2.51ogj'u + 
M Qj u + 2.51ogy'v — Ma y, where j\j and y'v are computed by adding the lu- 
minosities of the individual galaxies. The two methods return very similar 
values for the U — V colors. 



To convert between the measured rest-frame U — V color 
and the mass-to-light ratio M/L, we have generated stellar 
population synthesis models with the evolutionary synthesis 
code developed by G. Bruzual and S. Chariot (Bruzual & 
Chariot 2003). We selected the "Padova 1994" evolution- 
ary tracks, which are preferred by Bruzual & Chariot over 
the more recent "Padova 2000" tracks because the latter may 
be less reliable and predict a hotter red giant branch leading 
to worse agreement with observed galaxy colors. We used 
the solar metallicity set of tracks. The metallicities of the 
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TABLE 7 

DRGs/Red Galaxies Contribution to the Global Luminosity and Stellar Mass Densities 



Redshift 



Rest-frame Band 



Sample 



bf 
JP* 



2.0<z<2.5 


R 


J- 


-K>2.3 


a 07+0.04 
u - z/ -0.05 


n ?Q +<)03 

u - zy -0.03 


n4 ,+0.05 
u - 4I -0.05 


0.62^;| 


0.66^ 


a fi7 +0.20 






U- 


-V>0.25 


f) AQ+0.06 


n <;2 +<) " 4 

u -- ,/, -Q.Q4 


UDy -0.05 


0.84**> 


0.87|? 


g2+0.05 
u - yz -0.13 


2.7<z<3.3 


V 


J- 


-K>23 


A 7Q+0.06 

u - zy -0.05 


u - J:, -0.03 


a 44.+O.O6 
u - 4 ^-0.05 


u - Dy -0.28 




u - /y -0.22 






u- 


-V>0.25 


a 27 +0.04 
u - z/ -0.06 


u - J1 -0.04 


a 44.+O.O6 
u - 4 ^-0.05 


u - D;, -0.30 




79+*- 13 
u - /y -0.25 


2.0<z<2.5 


B 


J- 
u- 


-K>23 
-V>0.25 


1 9+o.os 
a , 7 +8:85 


u - zz -0.03 
u.tj_ 004 


a 31+0.06 

-Q-Q5 

A ro+0.06 

U - 5S -0.07 


62 m26 


a A7 +0.21 
U - D/ -0.31 

Q.89*fj 


0.74^; 2 ° 
0.90&» 


2.5<z<3.5 


B 


J- 


-K>2.3 


n 26+ 012 
u - zo -o.io 


u - JU -0.03 


a 30+0.05 
-»•!» 




0.82pj 


0.85$! 






u- 


-V>0.25 


a 2^+0.12 


29+0' 05 
U - Z '-O.03 


u '- ,o -0.06 


0.76*1 


o.si^if 


0-843:1 



en 
O 




NOTE. — The quoted errors correspond to the 1 a errors. 

a / PL is the contribution of DRGs/red galaxies to the global luminosity density; the three values correspond to the luminosity limits down 
to which the LF has been integrated: no limit, last observed point, and deep NIR MUSYC limit, respectively; see footnote of Table 6 for the 
specific integration limits. 

b f Pt is the contribution of DRGs/red galaxies to the global stellar mass density as estimated in § 5.3; the three values correspond to the 
luminosity limits down to which the LF has been integrated to estimate the luminosity density. 

1 Ci 

DRGs are poorly known, with evidence for solar and super- u 
solar metallicities for luminous DRGs (van Dokkum et al. 
2004). These DRGs appear more metal-rich than the five 
LBGs at z ~ 3 studied by Pettini et al. (2001, 0.1-0.5 Z ) 
and similar to the seven UV-selected star-forming "BX/MD" 
objects at z ~ 2 for which Shapley et al. (2004) inferred so- 
lar, and possibly supersolar, metallicities. In all cases, how- 
ever, the determinations rely on limited samples and suffer 
from large uncertainties. As shown by Ftirster Schreiber et 
al. (2004) adopting subsolar metallicity (Z = O.2Z ), the esti- 
mated M/L ratios are systematically lower by a factor of w 2 
on average. Therefore, if non-DRGs (or blue galaxies) are 
characterized by lower metallicities with respect to DRGs (red 
galaxies), the differences in M/L (and stellar mass densities) 
would be even larger than what is estimated assuming solar 
metallicities for both subsamples. For the star formation his- 
tory (SFH) we used three different prescriptions: a constant 
star formation history (CSF model), an exponentially declin- 
ing in time SFH characterized by the parameter r (tau-model), 
and an instantaneous burst model (SSP model). Several val- 
ues of r were used, from t = 0.1 Gyr (the resulting model 
being similar to the SSP model) to 6 Gyr (closer to the CSF 
model). We adopted the Chabrier (2003) IMF with lower 
and upper IMF mass cutoffs wl =0.1 M fi and m\j = 100 M , 
respectively. Adopting a different IMF would result in differ- 
ent derived mass-to-light ratios, which strongly depend on the 
shape and cutoff of the low-mass IMF (for example, assum- 
ing a Salpeter [1955] IMF, the estimated M/L ratio would be 
systematically larger by a factor of ~ 1 .7). However, since we 
are interested only in the relative contribution of DRGs (red 
galaxies) to the global stellar mass density, the results do not 
depend on the adopted IMF as long as all galaxies are char- 
acterized by the same IMF. We assumed that the interstellar 
extinction by dust within the objects followed the attenuation 
law of Calzetti et al. (2000) derived empirically from observa- 
tions of local UV-bright starburst galaxies under the formal- 
ism of a foreground screen of obscuring dust. 

We plot in Figure 9 the relation between the rest-frame 
U-V color and the mass-to-light ratio in the rest-frame R 
band, M/Lr, for the generated model tracks in the two cases 
with no extinction and with Ay = 1 . It is seen that dust extinc- 
tion moves the tracks roughly parallel to the model tracks. As 
emphasized by Bell & de Jong (2001), dust is a second-order 
effect for estimating stellar M/L ratios. Dust extinguishes 



1 
U-V 

FIG. 9. — Relation between the rest-frame U-V color and the mass-to-light 
ratio in the rest-frame R band, M/Lr, for the model tracks described in § 5.3. 
The dotted lines are for models with Ay = 0, while the other tracks are for 
Ay = 1- From top to bottom, the solid lines correspond to the CSF model 
(green), the tau-models (t = 6, 3, 1, 0.3, 0.1 Gyr) (red), and the SSP model 
(black). The blue dashed line, used in § 5.3 to estimate the contribution of 
DRGs and red galaxies to the total stellar mass density, corresponds to the 
median value of the tracks with A v = 1 ■ The arrow indicates the vector used 
to redden the Ay = model to Ay = 1 model. 

light from the stellar population, making it dimmer; however, 
dust also reddens the stellar population, making it appear to 
have a somewhat larger stellar M/L ratio. To first order, these 
effects cancel out, leaving a dust-reddened galaxy on the same 
color-stellar M/L ratio correlation. 

Using the relation between color and M/L, we convert the 
estimated median rest-frame U-V color and the measure- 
ments of the luminosity densities p L to stellar mass density es- 
timates p*. Specifically, we adopted the median value of M/L 
within the family of considered model tracks with Ay = 1 
(blue dashed line plotted in Figure 9 for the rest-frame R 
band); the error on the M/L was chosen as half of the dif- 
ference between the upper and lower envelope of the model 
tracks. The median rest-frame U-V color of each subsam- 
ple and the corresponding M/L in the rest-frame B, V, and R 
bands are listed in Table 8. 

Although the values of the stellar mass densities of the indi- 
vidual subsamples might be affected by very large uncertain- 
ties, the relative contribution to the global stellar mass density 
of DRGs (red galaxies) and non-DRGs (blue galaxies) should 
be more robust. The contribution f Pir of the DRGs (red galax- 
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TABLE 8 

Rest-Frame U-V Colors and EstimatedM/L Ratios 



Redshift 


Rest-Frame Band 


Sample 


a j/-y 


b logM/L 


Z.U Z.J 


n 
A 


J — K > Z.J 


U.OD ^U.O/ ) 


nn^J-n ~11 (f\ (\f\ -4- f) 






J — A \ Z.J 


— U.Uu (uAy) 


— U.oj±U.U/ l_-U.41zhU.lU; 






u — v ^ U.Zj 


U.OU (^J.OZJ 


n nn _i_ n. 77 en n 1 _i_ n 7 8^ 

— U.UU±U.Z/ (U.Ul±U.ZoJ 






U — V <_ U.Zj 


—U.ZJ (— U.1JJ 


fl 70 _|_ O flA (■ 71-1-0 OQA 
— U. /o ± U.U4 ^— U. / 1 ± \)X)y) 


2.7 < z < 3.3 




J—K > 2.3 


^7 CO ^7 ,> 








/-AT < 2.3 


-0.19 (-0.14) 


-0.80 ±0.05 (-0.75 ±0.05) 






{/-V >0.25 


0.59 (0.57) 


-0.05 ±0.30 (-0.08 ±0.28) 






{/-V <0.25 


-0.18 (-0.14) 


-0.75 ±0.05 (-0.75 ±0.05) 


2.0<z<2.5 


B 


7-A">2.3 


0.66 (0.67) 


0.00 ±0.40 (0.00 ±0.40) 






7-A"<2.3 


-0.06(0.19) 


-0.85 ±0.05 (-0.60 ±0.15) 






U-V >0.25 


0.60 (0.62) 


-0.10 ±0.35 (-0.05 ±0.35) 






{/-y <0.25 


-0.23 (-0.13) 


-1.15±0.15 (-0.85 ±0.05) 


2.5<z<3.5 


B 


7-AT>2.3 


0.56 (0.56) 


-0.13±0.33 (-0.13±0.33) 






7-a"<2.3 


-0.26 (-0.20) 


-1.15±0.15 (-1.10±0.20) 






C/-y> 0.25 


0.57 (0.57) 


-0.13±0.33 (-0.13±0.33) 






{/-y <0.25 


-0.25 (-0.20) 


-1.15±0.15 (-1.10±0.20) 



"Median rest-frame U-V color of the considered subsample; the numbers in parentheses correspond to the median U-V 
colors of the subsamples brighter than the MUSYC magnitude limit. 

b logM/L is the logarithm of the mass-to-light ratio in solar units assuming Ay = 1 mag, calculated as described in § 5.3; the 
numbers in parentheses correspond to the mass-to-light ratios of the samples brighter than the MUSYC magnitude limit. 



ies) to the global stellar mass density is listed in Table 7. 

Adopting the same assumptions for the stellar population 
synthesis models of DRGs (red galaxies) and non-DRGs (blue 
galaxies) (i.e., the median value of the considered model 
track), we see from Table 8 that DRGs (red galaxies) have 
M/L ratios systematically higher than non-DRGs (blue galax- 
ies) by a factor of ~ 5-11 depending on the rest-frame band 
(higher in the bluer bands). The differences in M/L are 
smaller when the brighter samples (down to the MUSYC 
limit) are considered, with the M/L of DRGs (red galaxies) 
being a factor of <~ 3-9 (4-9) larger with respect to non-DRGs 
(blue galaxies). For comparison, from the analysis of Spitzer- 
IRAC imaging on HDF-S, Labbe et al. (2005) found that the 
average mass-to-light ratio of DRGs in the rest-frame K band 
is about a factor of <~ 3 larger than LBGs, finding a corre- 
lation between M/L& and rest-frame U-V color. Consis- 
tently, from SED modeling of the 2 < z < 3.5 FIRES galax- 
ies, Forster Schreiber et al. (2004) found that the median rest- 
frame V-band M/L of DRGs is - 1 .2-2.3 M L^ (~ 0.4 for 
LBGs). The higher values of M/L for DRGs agree very well 
also with the results from SED fitting of individual galaxies 
in van Dokkum et al. (2004), who found for DRGs a median 
value M/Ly w 0.8 (M/L)©. Finally, our estimated M/L v are 
in excellent agreement with those estimated in Rudnick et al. 
(2006). 

Because of the systematically larger mass-to-light ratios, 
DRGs (red galaxies) dominate the global stellar mass density, 
with contributions in the range 66%-82% (69%-89%) down to 
the faintest probed rest-frame luminosities. The contribution 
of DRGs (red galaxies) increases up to 67%-85% (79%-92%) 
if the brightest sample is considered (down to the rest-frame 
magnitude limits of the deep NIR MUSYC). These numbers 
are consistent with the results of van Dokkum et al. (2006). 
From a complete mass-selected sample (M* > 10 11 M Q ) con- 
structed with MUSYC, FIRES, and GOODS-CDFS, they es- 
timated that DRGs in the redshift interval 2 < z < 3 make up 
77% of the total stellar mass, in very good agreement with 
our results. Our results are in qualitative good agreement also 
with the results from Rudnick et al. (2006), who found that 
the DRGs contribute 64% of the stellar mass density at z ~ 2.8 



and 30%-50% at z ~ 2. 

We stress that the estimated contributions of DRGs and red 
galaxies to the global stellar mass density are very uncertain 
and need confirmation from detailed SED analysis of mass- 
limited (rather than luminosity-limited) samples. 

6. SUMMARY AND CONCLUSIONS 

In this paper we have measured the rest-frame optical (R, 
V, and B band) luminosity functions of galaxies at redshifts 
2 < z < 3.5 from a composite sample constructed with the 
deep NIR MUSYC, the ultra-deep FIRES, and the GOODS- 
CDFS. The large surveyed area (<~ 378 arcmin 2 , 76% of 
which comes from the deep NIR MUSYC) of the compos- 
ite sample and the large range of luminosities spanned allows 
us to measure the bright end of the LF and to constrain the 
faint-end slope. Moreover, the several independent fields and 
their large area enabled us to largely reduce uncertainties due 
to sample variance, especially at the bright end. We have used 
Monte Carlo simulations to show that the uncertainties in the 
photometric redshift estimates do not significantly affect the 
measured parameters of the LF in the studied redshift regimes. 
There is a hint for a steepening in the faint-end slope of the LF 
from the rest-frame R band to the B band, although the differ- 
ences are not significant. The measured LF faint-end slopes at 
z > 2 are consistent, within the errors, with those in the local 
LFs. The characteristic magnitudes are significantly brighter 
than the local ones (e.g., ~ 1 .2 mag in the rest-frame R band), 
while the measured values for $* are typically a factor of ~ 5 
smaller with respect to the local values. 

The large number of objects in the composite sample al- 
lowed the first measurement of the LF of DRGs (defined 
based on their observed J-K color), which we compared to 
that of non-DRGs in the same redshift range. The DRG popu- 
lation is characterized by a very different LF than that of non- 
DRGs, especially at the faint end. While at the bright end 
the LF of DRGs is similar to that non-DRGs, at the faint end 
the latter one has a significantly steeper faint-end slope, espe- 
cially in the rest-frame R band. The significance of the dif- 
ference between the LFs of DRGs and non-DRGs decreases 
going to bluer rest-frame bands and to higher redshifts, al- 
though this is mainly caused by decreasing constraints on the 
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faint end of the LF of DRGs. Qualitatively similar results are 
found if we compare the LFs of red (rest-frame U— V > 0.25) 
and blue (rest-frame U — V < 0.25) galaxies in the same red- 
shift intervals, with the former equally contributing (or even 
dominating) at the bright end and the latter dominating the 
faint end. 

In the rest-frame V band we have also compared the LFs of 
blue galaxies (non-DRGs) with those of LBGs in the same 
redshift range. Although the two LFs agree very well at 
the bright end, the faint-end slope estimated by Shapley et 
al. (2001) is much steeper than the one measured in this pa- 
per. As the rest-frame optical LF of LBGs was estimated in 
Shapley et al. (2001) from the rest-frame UV LF and the ob- 
served distribution of &-K s colors as a function of 3?, mag- 
nitude, their steeper slope could be a result of an overesti- 
mate of the regression slope of the correlation between 
K s and 2% and/or of the faint-end slope of the adopted rest- 
frame UV LF of LBGs. Supporting the former possibility is 
the work of Labbe et al. (2006), who do not find any positive 
correlation between R magnitudes and R-K s colors in their 
deeper sample. Support for the latter comes from the very 
recent work of Sawicki & Thompson (2006), who measured 
the rest-frame UV LF for z ~ 3 LBGs from the Keck Deep 
Fields (KDF; Sawicki & Thompson 2005) and find a faint- 
end slope a = -1 .43^(j Qg, significantly shallower than the one 
adopted in Shapley et al. (2001). Alternatively, our blue K s - 
selected (i.e., rest-frame optical selected) galaxies might sim- 
ply constitute a different population than the z ~ 3 LBGs (rest- 
frame UV selected), with different characterizations of the LF 
at the faint end. We also caution that our measurements of 
the faint-end slopes still have significant uncertainties due to 
small number statistics. 

We generally find good agreement between our measured 
rest-frame B-band LFs at 2.5 < z < 3.5 and those previously 
published by Poli et al. (2003), Giallongo et al. (2005), and 
Gabasch et al. (2004). In the redshift range 2 < z < 2.5, the 
agreement between our rest-frame B- and /?-band LFs and 
those measured by Gabasch et al. (2004, 2006) from the FDF 
survey is less good, especially for the R band. Their Schechter 
parameters a and M* are consistent with ours only at the 2 er 
level in the rest-frame R band, while their estimated <£>* is 
larger than ours by a factor of ~ 1.3-1. 6 in the rest-frame B 
and R band, respectively. We have shown that this disagree- 
ment may be due to the spectroscopically confirmed existence 
of a (proto)cluster at z = 2.35 in the single field FDF survey. 

From the measured LFs we have estimated the number and 
luminosity densities of the global population of high-z galax- 
ies and of various subsamples. The contribution of DRGs (red 
galaxies) to the global number density is only ~ 13%-25% 
(18%-29%) down to the faintest probed rest-frame luminosi- 
ties. However, we have shown that field-to-field variations 
can be very significant (up to a factor of ~ 3), especially for 
relatively bright samples, in accord with the highly clustered 
nature of high luminosity and red galaxies (Adelberger et al. 
2005; Daddi et al. 2003; Quadri et al. 2007a). The contribu- 
tion of DRGs (red galaxies) to the global luminosity density 
is - 20%-30% (30%-50%), higher in the redder rest-frame 
bands (which are less affected by extinction and better tracers 
of the underlying stellar mass) and at lower redshifts. With 
respect to the lower z luminosity density estimates from the 
literature, we confirm the trend of slowly decreasing rest- 
frame /?-band luminosity densities beyond z <~ 1, with p L R 
at z ~ 2.3 being a factor of ~ 2 smaller than the local one. In 
the rest-frame B band, the measured global luminosity density 



at z ~ 2.3 is similar to the local value. At z ~ 3, the estimated 
global luminosity density may be a factor of <~ 2 higher, sim- 
ilar to the values around z ~ 1 . 

Finally, using stellar population synthesis models, we have 
derived the mass-to-light ratios of the considered subsamples 
by converting the estimated median rest-frame U -V color 
into M/L. In the rest-frame R and V bands, the mass-to-light 
ratios of DRGs (red galaxies) are a factor of <~ 5 larger than 
non-DRGs (blue galaxies), consistent with previous works. In 
the rest-frame B band the difference in M/L is higher, up to 
a factor of <~ 11. Using the estimated M/L, we have quanti- 
fied the contribution of DRGs and red galaxies to the global 
stellar mass density, finding that the total stellar mass budget 
is dominated by DRGs (red galaxies), whose contribution is 
of order <~ 60%-80% of the global value. We caution that our 
M/L ratios estimates are very rough and characterized by very 
large uncertainties and need confirmation from detailed SED 
analysis. 

The main limitation of this work is the small number statis- 
tics at the very faint end of the LF, which is probed only by the 
ultra-deep FIRES. The faint-end slopes of the DRG and red 
galaxy subsamples are especially uncertain. To make further 
progress in the determination of the LF of different galaxy 
populations at z > 2 and to better constrain the global LFs, it is 
crucial to better probe the faint end of the LFs. This can only 
be achieved with ultradeep NIR imaging with high-quality op- 
tical data over many spatially disjoint fields, in order to im- 
prove the statistics at the faint end and to mitigate the effect of 
field-to-field variations. Although we have shown that well- 
behaved photometric redshift errors do not affect significantly 
the measurement of the LF, the heavy reliance on photomet- 
ric redshifts is another limitation of this work, since "catas- 
trophic" failures and systematic errors could potentially affect 
the LF measurements. Obtaining large numbers of spectro- 
scopic redshifts for /f-selected high-z sources has proven diffi- 
cult and extremely time consuming. Even though the success 
rate for measuring spectroscopic redshift for bright galaxies is 
high with NIR spectroscopy (Kriek et al. 2006), only the use 
of multiobject NIR spectrographs will make it possible to con- 
struct a large sample of high-z ^-selected galaxies with spec- 
troscopic redshift measurements. Further advances can be ex- 
pected from further extension of the wavelength range into 
the red. Scheduled Spitzer IRAC observations on the deep 
NIR MUSYC fields will allow us to (1) separate old and pas- 
sively evolving galaxies from heavily obscured and strongly 
active star-forming galaxies (see Labbe et al. 2005), making it 
possible to study the LF of physically different types of galax- 
ies; (2) extend the robust measurement of the LF at redshift 
2 < z < 3 into the rest-frame NIR, which is much closer to a 
selection by stellar mass; (3) convert the measured rest-frame 
NIR luminosity function into a mass function and study the 
evolution of the stellar mass density; and (4) extend the study 
of the rest-frame optical LFs to even higher redshifts. 
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FIG. Al . — Left: Intrinsic Schechter LF (red dashed line) compared to the observed LF at 0.5 < z < 0.8 corresponding to the median Monte Carlo realization 
(black solid line and filled circles correspond to the maximum likelihood and 1 /V max method, respectively); the best-fit M* are also plotted as arrows. Right: 
Measured Schechter parameters M* plotted vs. a; the black filled circle represents the input parameters of the "observed" model catalog, together with the 1, 2, 
and 3 a contour levels; the green filled triangles represent the 100 Monte Carlo realizations; the red filled square is the median realization, and the plotted error 
bars represent the rms in a and M* . The amount of the systematic effects in a and M* are also specified. The red filled star represents the median realization 
when the redshifts of the objects are extracted from a random uniform distribution, as done in Chen et al. (2003); the small black filled circles represent the 
corresponding 100 Monte Carlo realizations. The result from Chen et al. (2003) is also shown (blue open star) together with its 99% uncertainties (blue dashed 
contour). 
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APPENDIX A 



EFFECTS OF UNCERTAINTIES IN Z PHO t: MONTE CARLO SIMULATIONS 

In order to quantify the systematic effects on the LF parameters a and M* due to the uncertainties in the photometric redshift 
estimates, we performed a series of Monte Carlo simulations. First, we generated several model catalogs of 25,000 galaxies with 
redshifts between Zi mc an d Z2,mc an d with luminosities drawn from an input Schechter LF. While in the Monte Carlo simulations 
of Chen et al. (2003) the redshifts of the objects in the mock catalogs were extracted from a random uniform distribution, we took 
into account the fact that, under the assumption of no evolution in the number density, the probability of a galaxy existing at the 
redshift z is proportional to the volume: 

P M - dv - d M 1 (A1) 

P( ) dz (1 + z) 2 ^a + «m(1+z) 3 ' 

where dt is the luminosity distance. Since in a real survey galaxies are selected down to a limiting apparent magnitude, the final 
mock catalogs are obtained after applying a limit in the observed apparent magnitude. The effect of a limiting apparent magnitude 
is that, at a fixed observed magnitude, intrinsically fainter sources are systematically excluded from the catalog at higher redshift. 
Next, as done by Chen et al. (2003), we assumed a redshift error function parameterized as a Gaussian distribution function of 1 
a width tj' z (\+z), with a' z the scatter in Az/(1 +z spe c), and we formed an observed redshift catalog by perturbing the input galaxy 
redshift within the redshift error function. Finally, we determined the LF for the galaxies at zi < z < Z2 using the 1 /Vmax and 
maximum likelihood methods described in § 3. Note that we ignore /^-correction in our Monte Carlo simulations. 

We first studied the effects of the photometric redshift uncertainties at z < 1, by using zi.mc =0.15, Z2,mc = 1-6, and assuming 
an input Schechter LF with parameters a = -1 .2 and M* = —21 .5 (as in Chen et al. 2003). In order to compare our results with 
those from the Monte Carlo simulations in Chen et al. (2003), we used a' z = 0.15 and we measured the LF in the same redshift 
range 0.5 < z < 0.8. We find that the median measured M* is brighter than the intrinsic value by ~ 0.4 mag and that the measured 
a is steeper on average than the intrinsic value by ~ 0.1. This result is shown in Figure Al: in the left panel, the input Schechter 
LF is compared to the median Monte Carlo realization; in the right panel, the 100 Monte Carlo realizations are plotted in the 
a-M* plane and compared to the best-fit values and the corresponding 1, 2, and 3 a contour levels of the LF measured on the 
redshift-unperturbed mock catalog. 

As shown in Figure Al, the measured systematic effects on a and M* caused by the redshift uncertainties arise from an excess 
of sources at both the faint and the bright end of the LF. A careful analysis of the Monte Carlo simulations reveals the origin of 
these excesses. Because of the uncertainties in the redshifts, sources can scatter from high to low redshifts, and vice versa. As 
the probability of a galaxy existing at redshift z is proportional to dV /dz (which peaks at z ~ 2.5), the number of sources that 
scatter from higher redshifts to lower ones is much larger than vice versa. This is evident in Figure A2, where the rest-frame 
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FIG. A2. — Rest-frame absolute magnitudes M plotted vs. the redshifts z of the redshift-unperturbed mock catalog (gray points). The solid line represents the 
limiting magnitude of the sample as a function of redshift; the two vertical dashed lines represent the redshift interval within which the LF is measured. Cyan 
filled circles at z > 0.8 are the sources that in the redshift-perturbed mock catalog scatter from high z into the considered redshift bin (red filled circles); the cyan 
filled circles at z < 0.5 are the sources that scatter from low z into the considered redshift bin (green filled circles). The loci at magnitudes fainter than the lower 
dotted line and brighter than the upper dotted line are the regions where we observe the excess of sources at the faint and bright ends (evident in the left panel of 
Fig. Al), respectively. 

absolute magnitudes of the mock catalog are plotted versus redshift: at a fixed magnitude, the number of sources is larger at 
higher redshifts. 

Since an object in the mock catalog is characterized by a fixed apparent magnitude, a new redshift estimate translates into a 
new rest-frame absolute magnitude; e.g., when a source "scatters" from high to low redshift, the estimated absolute brightness 
is fainter than the intrinsic value, and the object moves, in Figure A2, from right to left along a line parallel to the black solid 
line. Moreover, at any redshift, there are more sources at faint magnitudes than at bright ones because of the shape of the LF. 
Therefore, the sources scattering from high redshifts into the considered redshift bin (plotted in Fig. A2 in cyan) preferentially 
end up at fainter magnitudes (plotted in red), producing the excess at the faint end of the LF with respect to the input LF. The 
excess at the bright end is instead mainly produced by those sources that scatter into the considered redshift range coming from 
lower redshifts (represented in Fig. A2 with cyan and green filled circles); since at the bright end the number of objects is very 
small, even a handful of new sources can significantly increase the measured density with respect to the intrinsic one. Because at 
low redshifts the dependency of dM/dz with redshift is strong, even small uncertainties in the redshift estimate have large effects 
on the rest-frame absolute magnitude. For example, at z ~ 0.7, a Az — 0.25 results in a AM ~ 0.7- 1 mag (only 0.1-0.2 mag at 
z = 2.25). 

The systematic effect on M* that we derive (AM* ~ -0.4 mag) is about half the effect found by Chen et al. (2003). Also 
almost no systematic effect in a was found in their work. It seems very likely that these differences might be due to the fact that 
the redshifts are drawn from a random uniform distribution in Chen et al. (2003) while in our Monte Carlo simulations they are 
extracted from the probability function specified in Eq. Al. In the former case, there would be a much larger number of low-z 
sources that can scatter into the considered redshift bin, resulting in a larger excess of sources at the bright end and therefore a 
larger systematic effect in M*. If we repeat our Monte Carlo simulations extracting the redshifts of the sources from a random 
uniform distribution, we obtain AM* ~ -0.9 mag and Act ~ -0.06, consistent with the result of Chen et al. (2003) as shown in 
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FIG. A3. — Left: Input Schechter LF (red dashed line) compared to the observed LF at 2 < z < 2.5 of the median Monte Carlo realization (black solid line and 
filled circles correspond to the maximum likelihood and 1 /V max methods, respectively); the best-fit M* are also plotted as arrows. Right: Measured Schechter 
parameters M* plotted vs. a; the filled black circle represents the input parameters of the "observed" model catalog, together with the 1, 2, and 3 cr contour levels; 
small green filled circles represent the 100 Monte Carlo realizations using a' L = 0.12; the red filled star is the corresponding median realization; the red filled 
triangle and square represent the median realization of the Monte Carlo simulations with a' z = 0.05 and a' z = 0.3, respectively. The amount of the systematic 
effect in a and M* is also specified for cr' z = 0.12. 



the right panel of Figure Al . 

Next, we repeated our Monte Carlo simulations at higher redshift by generating model catalogs with galaxies at redshifts 
between z^mc = 1 an d ^2.mc = 6, since the goal of this paper is to measure the LF of galaxies in the redshift intervals 2 < z < 2.5 
and 2.5 < z < 3.5. In Figure A3 we plot the results of our Monte Carlo simulations at 2 < z < 2.5, assuming an input Schechter LF 
with a = -1 .3 and M* = -22 and for a' z =0.12 (which corresponds to the photometric redshift errors in the deep NIR MUSYC 
forz > 1.5 objects). 

The systematic effects on the measured a and M* are now very small, AM* = -0.05 mag and Act = -0.007, and negligible with 
respect to the other uncertainties on the estimated best-fit parameters. Similar results are obtained in the redshift bin 2.5 < z < 3.5. 
At 2 < z < 3.5, the effect due to p(z) °= dV jdz is much smaller than at z < 1, since p(z) peaks at z ~ 2.5 and then decreases very 
slowly, so that the number of high- and low-z sources scattering into the considered redshift bin is similar to the number of sources 
scattering out. Also, at z > 2, the error on the rest-frame absolute magnitude corresponding to a redshift error is significantly 
smaller than at z < 1; e.g., AM = 0.3-0.4 mag for Az ~ 0.4 at z — 2.25. Therefore, the measured LF is similar to the input one 
and the systematic effects on a and M* are negligible compared to the uncertainties in the LF estimates for reasonable values of 
OS 0.12). 

We repeated the above Monte Carlo simulations assuming different input a (-1.5, -1, and -0.5) to study the behavior of 
the systematic effects as function of the faint-end slope. No significant differences are found: -0.04 < AM* < 0.01 mag for 
a = -1 and -0.5 (depending on the considered redshift bin); for a = -1.5, the systematic effect is slightly larger (AM* ~ 
-0.07/ + 0.08 mag for 2 < z < 2.5 and 2.5 < z < 3.5, respectively), but also the uncertainties on the best-fit M* increase with 
steeper faint-end slopes (since the observed LF appears more like a power law) so that the systematic effects on the measured 
best-fit Schechter parameters remain very small with respect to the uncertainties on the best-fit values. 

Finally, we investigated the effects of non-Gaussian redshift error probability distributions. First, using a model catalog with 
galaxies at redshifts between zi.mc =0.1 and Z2,mc = 6, we simulated the effect of a 5% "catastrophic" outliers by assigning 
random redshifts to 5% of the mock catalog. Adopting the input LF with a = -1 .3 and M* = -22.0 and assuming a' z =0.12, we 
find larger systematic effect in both a and M* by a factor of almost 2. Next, we built the mean redshift probability distribution 
of the deep NIR MUSYC by averaging the individual redshift probability distribution for each galaxy calculated by the used 
photometric redshift code (for details see Rudnick et al. 2003). The average MUSYC redshift probability distribution is well 
modeled by a Lorentzian function, rather than a Gaussian function. We find a larger systematic effect in a, twice as much as the 
corresponding effect assuming a Gaussian parametrization for the redshift probability distribution, but similar systematic effect 
in M*, although in the opposite direction. To summarize, although the systematic effects in a and M* expectedly get larger when 
we simulate "catastrophic" outliers or we adopt a redshift error function with broader wings compared to the Gaussian model, 
they remain much smaller than the random uncertainties in the LF estimates. 

We also quantified the systematic effect on the luminosity density estimates. We find that the effect is of the order of a few 
percent (always <6%) depending on the considered redshift interval and on the input Schechter LF. 

In order to include this contribution in the error budget, we conservatively assume a 10% error contribution to the luminosity 
density error budget due to uncertainties in the photometric redshift estimates. 
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FIG. Bl.— Rest-frame V-band LF at 2.7 < z< 3.3. Top: Rest-frame V-bandLFs at 2.7 < z< 3.3 for DRGs/non-DRGs (left) and red/blue galaxies (right). The 
red solid line and filled circles represent the LF of DRGs and red galaxies estimated with the STY and 1 /V max methods, respectively; the blue dot-dashed line and 
filled squares represent the LF of non-DRGs and blue galaxies. The rest-frame V band LF of all galaxies is also plotted (black dashed line). The shaded regions 
represent the 1 a uncertainties of the LFs measured with the STY method. The arrows represent the best-fit M*\ error bars as in Fig. 5. Bottom: 1, 2, and 3 a 
contour levels from the STY method for DRGs and red galaxies (red dashed lines), for non-DRGs and blue galaxies (blue solid lines), and for all galaxies (black 
dot-dashed lines); the best-fit Schechter values are also plotted (red filled circle for DRGs and red galaxies, blue filled square for non-DRGs and blue galaxies, 
and open circle for all galaxies). 

APPENDIX B 

REST-FRAME V- AND fi-BAND LUMINOSITY FUNCTIONS FOR DRGS/NON-DRGS AND 

RED/BLUE GALAXIES 

In the body of the paper we showed the rest-frame R band LFs of DRGs, non-DRGs, red and blue galaxies (see Fig. 5). For 
completeness, we show here the comparison of the LFs of DRGs (red galaxies) and non-DRGs (blue galaxies) discussed in § 4.2 
in the rest-frame V band at 2.7 < z < 3.3 (Fig. Bl) and in the rest-frame B band at 2 < z < 2.5 (Fig. B2) and at 2.5 < z < 3.5 
(Fig. B3). The corresponding best-fit Schechter parameters are listed in Table 3. 

APPENDIX C 

COMPARISON WITH PREVIOUSLY PUBLISHED LUMINOSITY FUNCTIONS 



Here we compare our results to previous rest-frame optical LF studies at z > 2, which were based on smaller samples. We note 
that these studies are affected by significant uncertainties due to field-to-field variance (as they are based on a single field or on a 
very small total surveyed area) and by small number statistics at the bright end. 




CI. POLIETAL. (2003) AND GIALLONGO ET AL. (2005) 

Poli et al. (2003) analyzed a sample of 138 /f-selected galaxies down to A'ab = 25 to construct the rest-frame B-band LF in 
the redshift range 1.3 < z < 3.5. The total area of their composite sample is ~ 68 arcmin 2 , a factor of ~ 5.6 smaller than the 
area sampled in this work. Giallongo et al. (2005) repeated the analysis in Poli et al. (2003) with an improved composite sample 
(although with the same area) and allowing the Schechter parameters and M* to vary with the redshift, while a is kept constant 
at the low-redshift value (z < 1). A direct comparison between the LFs of Poli et al. (2003) and Giallongo et al. (2005) and the 
LF measured in this work is shown in Figure CI (left panel) for the redshift range 2.5 < z < 3.5. 

First, we note in Figure CI how much better the bright end of the LF is constrained from our work: the large area of the com- 
posite sample (~ 76% of which comes from the deep NIR MUSYC alone) allows us to sample the LF up to Z?-band magnitudes 
~ 1.5 mag brighter than done in Poli et al. (2003) and Giallongo et al. (2005). Our measurements of the LF using the 1 /y max 
method are consistent within the errors with those in both Poli et al. (2003) and Giallongo et al. (2005). The best-fit LF estimated 
with the maximum likelihood analysis in Poli et al. (2003) is consistent, within the errors, with our best-fit solution. However, 
the best-fit LF estimated with the maximum likelihood analysis in Giallongo et al. (2005) is significantly different from ours, as 
clearly shown in the inset of Figure CI (left panel). While their faint-end slope is consistent with our best-fit a, their lack of 
constraints on the bright end of the LF results in a much fainter Mg (by ~ 0.6 mag). 

C2. GABASCHETAL. (2004, 2006) 

Gabasch et al. (2004, 2006) analyzed a sample of 5558 /-selected galaxies down to Iab = 26.8 (50% compleness limit) from 
the FDF survey (Heidt et al. 2003) to study the evolution of the rest-frame B- and r-band LF over the redshift range 0.5 < z < 5. 




The total area of their sample is ~ 40 arcmin 2 over a single field, a factor of ~ 9.5 smaller than the total area sampled in this 
work. 

A direct comparison between the B-band LF from Gabasch et al. (2004) and ours is possible in the redshift range 2.5 < z < 3.5, 
and it is shown in the left panel of Figure CI (red symbols). The best-fit Schechter parameters a, M*, and <£>* are consistent with 
ours within the errors. In the redshift range 2 < z < 2.5 we compared our estimated Z?-band LF with the one defined by their 
best-fit Schechter parameters (a, M*, and $*) estimated at z ~ 2.25. While their best-fit Schechter parameters a and M* are 
consistent with ours, their best-fit $* is a factor of ~ 1.3 larger. This difference can be entirely accounted for by field-to-field 
variations (see §5.1). In fact, the FDF survey consists of a single pointing of only ~ 40 arcmin 2 and thus it is potentially strongly 
affected by sample variance. Moreover, Gabasch et al. (2004) spectroscopically identified an overdensity of galaxies at z = 2.35 
(possibly a [proto]cluster, with more than 10 identical redshifts), which can potentially strongly bias the estimate of $* in this 
redshift bin. 

In the right panel of Figure CI we compare our measured /?-band LF with that measured by Gabasch et al. (2006) in the 
redshift range 1.91 <z<2.61. The agreement between the two LFs is much worse than for the rest-frame fi-band. Their best-fit 
Schechter parameters a and M* are now consistent with ours only at the ~ 2 a level, with their M* about 0.3 mag brighter 
than ours; their best-fit $* is a factor of ~ 1.6 larger than our best-fit value. As mentioned above, an overdensity of galaxies 
was spectroscopically found in this redshift interval in the FDF. Although it is hard to quantify the effect of the presence of a 
(proto)cluster at this redshift on the measured LF, it is interesting to make a connection with the work of Steidel et al. (2005), who 
spectroscopically identified a protocluster at z = 2.3 in the HS 1700+643 field. In the spectroscopically identified galaxy sample, 
Steidel et al. (2005) found 19 (out of 55) objects at z = 2.300 ±0.015 in the redshift range 2 < z < 2.5. Within the - 7' x 7' field 
of view (similar to the FDF field of view) over which the protocluster objects are distributed, the protocluster and "field" galaxy 
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Fig. Cl . — Left: Rest-frame B-band LF at 2.5 < z < 3.5 from this work (blue solid line and black filled circles) compared to the LF from Poli et al. (2003; 
dotted line and green open squares), from Giallongo et al. (2005; dashed line and orange open circles), and from Gabasch et al. (2004; dot-dashed line and red 
open triangles). In the inset, the 1, 2, and 3 <r contour levels of the best-fitSchechter parameters a and M* are compared to the best-fit values estimated by Poli et 
al. (2003), Giallongo et al. (2005), and Gabasch et al. (2004). Right: Rest-frame S-band LF at 2 < z < 2.5 from this work (red solid line and black filled circles) 
compared to the LF from Gabasch et al. (2006; orange dashed line and filled triangles; R = r-0.12). In the inset, the 1, 2, and 3 o contour levels of the best-fit 
Schechter parameters a and M* are compared to the best-fit values estimated by Gabasch et al. (2006). 

sky distributions are the same. Assuming that the distribution of the spectroscopically identified galaxies is representative of the 
whole sample, the estimated density would be a factor of ~1. 3-1.6 that of "field" galaxies only, consistent with the difference 
found in $* between our composite sample and the FDF in the redshift range 2 < z < 2.5. 



C3. COMPARISON WITH GIALLONGO ET AL. (2005) FOR RED/BLUE GALAXIES 

Giallongo et al. (2005) measured the rest-frame fi-band LF of red and blue galaxies at redshift 2.5 < z < 3.5 using a sample of 
138 ^-selected galaxies in the redshift range 1.3 < z < 3.5 down to K^g = 25. Their red and blue populations were defined on 
the basis of an "SO color track". From their Figure 1, the average rest-frame U—V color of their model at redshift 2.5 < z < 3.5 
is ~ 0.22, very similar to our definition of blue and red galaxies (U — V < 0.25 and U — V > 0.25, respectively). In Figure C2 we 
have compared the LFs of red and blue galaxies from this work to those presented in Giallongo et al. (2005). 

As for the LF of all galaxies, we are able to constrain the bright end of the LF much better by sampling the LF to luminosities 
~ 1.5 mag brighter. For the red galaxy population, their measurements of the LF with the 1 /V ma x method are consistent within 
the errors with ours; their best-fit faint-end slope is also consistent, within the errors, with our estimate, although a was fixed at 
the local value in Giallongo et al. (2005); our M* is ~ 0.4 mag brighter than their best-fit value, but still consistent within the 
errors at the 1 a level. The best-fit normalization of Giallongo et al. (2005) is a factor of ~ 1.5 smaller than ours, only marginally 
consistent at the 2 a level. 

For the blue population, the 1 /V max measurements of the LF from Giallongo et al. (2005) are consistent with ours. However, the 
agreement gets worse when we compare our LF estimated using the STY method with theirs, especially at the bright end, where 
they lack information. While their assumed faint-end slope is consistent with our best-fit solution, the characteristic magnitude 
from Giallongo et al. (2005) is ~ 0.7 mag fainter than ours. Their best-fit solution is, however, very uncertain, characterized by 
very large error bars and still consistent with ours at the ~ 1 a level. 
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FIG. C2. — Left: Rest-frame B-band LF at 2.5 < z < 3.5 from this work for red galaxies (U — V > 0.25; red solid line and black filled circles) compared to the 
one from Giallongo et al. (2005; dashed line and orange open circles); in the inset, the 1, 2, and 3 a contour levels of the best-fit Schechter parameters a and 
M* are compared to the best-fit values estimated by Giallongo et al. (2005; orange open circle), which assumed a constant a = —0.46. Right: Rest-frame B-band 
LF at 2.5 < z < 3.5 from this work for blue galaxies (U — V < 0.25; blue solid line and black filled circles) compared to the one from Giallongo et al. (2005; 
dashed line and cyan open circles); in the inset, the 1, 2, and 3 a contour levels of the best-fit Schechter parameters a and M* are compared to the best-fit values 
estimated by Giallongo et al. (2005; cyan open circle), which assumed a constant a = —1.38. 
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